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A MONTE CARLO ERROR ANALYSIS PROGRAM FOR NEAR -MARS, 

FINITE -BURN, ORBITAL TRANSFER MANEUVERS 

By Richard N. Green, Lawrence H. Hoffman, 
and George R. Young 
Langley Research Center 

SUMMARY 

A computer program has been developed which performs an error analysis of a 
minimum -fuel, finite -thrust, transfer maneuver between two Keplerian orbits in the vicin- 
ity of Mars. The method of analysis is the Monte Carlo approach where each off-nominal 
initial orbit is targeted to the desired final orbit. The errors in the initial orbit are 
described by two covariance matrices of state deviations and tracking errors. The func- 
tion of the program is to relate these errors to the resulting errors in the final orbit. 

The equations of motion for the transfer trajectory are those of a spacecraft maneu- 
vering with constant thrust and mass -flow rate in the neighborhood of a single body. The 
thrust vector is allowed to rotate in a plane with a constant pitch rate. The transfer 
trajectory is characterized by six control parameters and the final orbit is defined, or 
partially defined, by the desired target parameters. 

The program is applicable to the deboost maneuver (hyperbola to ellipse), orbital 
trim maneuver (ellipse to ellipse), fly-by maneuver (hyperbola to hyperbola), escape 
maneuvers (ellipse to hyperbola), and deorbit maneuver. 

INTRODUCTION 

The effort of the United States to land an unmanned capsule on Mars is a very com- 
plex mission. The mission profile of Project Viking contains two or three midcourse 
maneuvers en route, a Mars orbit insertion maneuver, numerous orbital trim maneuvers 
for planetary photoreconnaissance and positioning of the spacecraft over the landing site, 
and finally a deorbit maneuver. The navigation and guidance problems are complicated 
further by the launch of a second spacecraft within a 50-day period. 

For a mission of this complexity it is necessary to determine the sensitivity of the 
trajectory with respect to error sources which requires the establishment of the expected 
magnitude of the various maneuvers. 



This paper deals with the error analysis of the near -Mars maneuvers (those within 
the sphere of influence of Mars). Of primary importance is the analysis of the Mars 
orbit insertion maneuver. This analysis is necessary in order to determine the orbit 
errors after the insertion maneuver, to calculate velocity budgets, and to determine the 
range of control variables for design purposes. 

Several assumptions about planetary orbital insertion maneuvers have been adopted 
in the past to simplify the error analysis. They are as follows: (a) the maneuvers can 
be analyzed impulsively, (b) the encounter errors can be mapped through the maneuvers 
using linear propagation theory, and (c) the encounter errors can be sufficiently defined 
by using the impact plane and time -of -flight parameters. The validity of these assump- 
tions is questionable for the Viking orbital insertion maneuver due to the long burn time 
required and the large approach errors caused by a long interplanetary trip time. There- 
fore, a computer program, VEAMCOP (Viking Error Analysis Monte Carlo Program), 
was developed to remove these assumptions from the analysis. The basic targeter for 
VEAMCOP is VITAP (Viking Targeting Analysis Program) given in reference 1. Although 
the mathematical structure of the program is quite general and therefore applicable to 
maneuvers near any planet, VEAMCOP was developed specifically for near -Mars maneu- 
vers and mainly for the orbit insertion maneuver. The maneuvers considered in the anal- 
ysis are minimum -fuel, finite-burn transfers between two conics. The guidance law 
allows for the targeting of the spacecraft to a final conic that is specified by certain con- 
straints. In addition, the law includes an optimum maneuver in the sense that the fuel 
required for the transfer is a minimum. 

The function of the program is to relate the errors in the initial orbit and errors 
in the thrusting maneuver to the resulting errors in the final orbit. The statistics of the 
errors in the final orbit are found by use of the Monte Carlo Method. First a sample 
estimate of the initial orbit is randomly generated from the statistics of the initial orbit. 
The random estimate is then targeted to the desired final orbit. However, due to the 
lack of knowledge of the actual initial orbit and due to maneuver execution errors, the 
final orbit will be in error. The Monte Carlo Method of repeatedly sampling the initial 
errors and calculating the final errors provides the statistics of the errors in the final 
orbit. 


SYMBOLS 

a semimajor axis, kilometers 

a(t) magnitude of thrust acceleration at time t, kilometers/second^ 
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B 

B-R 

B-T 

C 

D 

E(£) 

e 

F 

G 

i 

K 

M 

m 

xn(t) 

N s 

N(c,d) 

N 

N,V,W 


vector from the center of the planet perpendicular to the incoming hyper- 
bolic asymptote, kilometers 

targeting parameter (component of B in the R direction), kilometers 
targeting parameter (component of B in the T direction), kilometers 
covariance matrix 

Cartesian covariance matrix of actual state deviations from the nominal 
state 

expected value of 4 
eccentricity 

point which defines the computed initial direction of thrust (see fig. 6(a)) 
point which defines the actual initial direction of thrust (see fig. 6(a)) 
inclination, degrees 

number of velocity counting accelerometer pulses before thrust is 
terminated 

dimension of vector 

mass, kilograms 

mass of spacecraft at time t, kilograms 
number of samples 

normal density function with mean c and variance d 
unit vector in direction of ascending node 

spacecraft coordinate system where V points in direction of velocity, 

W is perpendicular to the orbital plane in direction of angular momentum 
vector, and N completes right-handed system 
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n unit vector normal to plane of thrust 

P,Q,W coordinate system (see fig. 6(b)) 

AAA 

p,q,w coordinates in P,Q,W-system 

R random number from normal distribution with mean 0 and variance 1 

R = S X T 

r radius from center of planet, kilometers 

r a radius of apoapsis, kilometers 

rp radius of periapsis, kilometers 

S unit vector parallel to incoming hyperbolic asymptote 

T Cartesian covariance matrix of tracking errors 

T unit vector in planet equator perpendicular to S 

t time, seconds 

t b time duration of thrusting maneuver, seconds 

U(f,g) uniform density function defined between f and g 

u random number from uniform distribution defined on interval [0.1] 

V integral of acceleration due to thrust 

V a actual velocity gained, kilometers/second (see eq. (6)) 

V b velocity counting accelerometer bias, meters/second 

Vtbg velocity -to-be-gained, kilometers/second (see eq. (5)) 

Voc hyperbolic excess velocity, kilometers/second 
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v velocity increment equivalent to one pulse of velocity count ing 

accelerometer, meters/second 

v cal calibrated value of v, meters/second 

AAA 

X,Y,Z rectangular Cartesian base vectors 

actual Cartesian state of spacecraft 
= X n + AX d R 
X e estimate of X a 

X n nominal Cartesian state of spacecraft 

x, y,z rectangular Cartesian coordinates, kilometers 

q!,/ 3,6 angles defining initial direction and plane of thrust, degrees (see fig. 1) 

y, p,ifs random variables, degrees (see fig. 6(a)) 

AX d = X a - X n 

AX e = X e - X n 
AX t = X e - X a 

e error in the calibration of the velocity counting accelerometer, meters/second 

t P 

rj random vector in principal axis system 

0 rotation matrix 

6,<fi angles of rotation, degrees (see fig. 6(b)) 

6 thrusting pitch rate, degrees/second (see fig. 1) 

X eigenvalues 

ju gravitational constant of planet, kilometers^/second^ 
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V 


true anomaly, degrees 


v Q true anomaly at start of maneuver, degrees 

£ arbitrary variable 

2 azimuth angle, degrees 

a standard deviation of random variable (symbols of variables will be used 

as subscripts to refer to the deviations of specific parameters) 

r magnitude of thrust, kN 

$ matrix of eigenvectors 

longitude of ascending node, degrees 

cu argument of periapsis, degrees 

Subscripts: 

D D -matrix 

T T -matrix 

i,j,k indices 

f final 

n nominal 

o initial condition 

x,y,z rectangular Cartesian coordinates 

1,2,. . . first, second, etc. 

Superscripts: 

j ,k indices 
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vector 


R random 

— R random vector 

— mean value 

T matrix transpose 

~ unit vector 

. differentiation with respect to time 

* specific value 

' ~y modified parameters 

Notation: 


~ distributed as 

[| £ || denotes integer part of | 

P r [x] probability of X occurring 

A increment 


ANALYSIS 

The problem considered in this paper is that of relating the statistics of the errors 
in the initial orbit and errors in the maneuver to the statistics of the errors in the final 
orbit where the final orbit results from a minimum fuel, finite burn, orbital transfer. 
First, the physics of the problem are formulated resulting in a set of three second-order 
differential equations which represent the equations of motion of the thrusting maneuver. 
Next, the targeting problem is considered which leads to the numerical values of the six 
control parameters that characterize the transfer from the initial orbit to the desired 
final orbit. The statistics of the initial orbit and the maneuver execution errors are then 
discussed. The Monte Carlo Method of analysis follows which consists of generating a 
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random sample of the initial orbit, targeting this sample to the desired final orbit, and 
computing the perturbed final orbit which results from errors in the initial orbit and 
execution errors. This sequence is repeated many times and the statistics of the errors 
in the final orbit are formulated. 


Equations of Motion 

The equations of motion for the transfer trajectory are those of a spacecraft maneu- 
vering with constant thrust r and mass -flow rate m in the neighborhood of a single 
body. If the mass of the spacecraft at the start of the maneuver is m 0 , then the mass 
at time t is given by 

m(t) = m Q + mt 

and the magnitude of the acceleration due to thrust by 

a(t) = ■“ t~t 

m 0 + mt 

The thrust vector is allowed to rotate in a plane with a constant pitch rate 6 so that the 
transfer trajectory is characterized by six control parameters a, /3, 6, 6, t^, v Q 
(see fig. 1). The orientation of the plane of thrust relative to an inertial X,Y,Z -coordinate 


Z 



Figure 1.- Sketch of coordinate system and angles used in equations of motion. 
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system is defined by the two angles a and 0. At the start of the maneuver the thrust 
vector is at an angle 6 from the X,Y plane and rotates in the a,i 3 plane at a rate 6 
until the burn is terminated at t = t^. The direction cosines of the thrust vector at time 
t are therefore given by 

cos(6 + 6t)cos a - sin(6 + 0t)sin a cos 0 
cos(6+ 0t)sin a + sin(5 + 0t)cos a cos 0 


sin(6 + 0t)sin 0 

The assumed trajectory model is two -body motion plus an acceleration due to thrust and 
is defined by the equations of motion 

x = - ^ + a(t)fcos(5 + 0t)cos a - sin(5 + 0t)sin a cos 0 ] 
r 3 -> 

y = -i^ + a(t ) cos(S + 0t)sin a + sin(5 + 0t)cos a cos 0 
r 3 L J 

z = - + a(t)[sin(6 + 6t)sin 0j 

where 


a(t) = 


m Q + mt 


= (x 3 + y^ + z 2 ) 1 ^ 2 


The initial conditions for the equations of motion are derived from the knowledge of the 
initial orbit plus the control parameter v Q which denotes the true anomaly on the orbit 
at the start of the maneuver. To conserve computer time the equations of motion are 
"integrated" by a truncated power series expansion in time (see appendix of ref. 1). 

Targeting Problem 

The guidance law requires that the control parameters be determined such that the 
fuel is minimized subject to the requirement that the final orbit satisfies certain con- 
straints. These constraints might consist of requiring the final orbit to have a specific 
orientation, period, and so forth. Actually, the constraints are chosen from a set of 
20 possible target parameters (ref. 1). Since the maneuver is defined by up to six con- 
trol parameters, at most six constraints can be satisfied. For the fuel to be minimized, 
the number of constraints must be less than the number of control parameters. The tar- 
geting problem, then, is to determine the set of control parameters which satisfy the 
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constraints imposed on the final orbit and minimize the fuel required for the transfer 
maneuver. Basically, the targeting involves the solution of a constrained minimization 
problem by use of constant Lagrange multipliers and the Newton -Raphs on iteration tech- 
nique. The targeter for VEAMCOP is VITAP (ref. 1). 

Statistics of the Problem 

The statistics of the errors in the initial orbit are defined by two covariance 
matrices. In the absence of all errors the spacecraft would proceed along the nominal 
trajectory. However, due to previous perturbations and error sources, the actual path 
of the spacecraft will not be the nominal trajectory but a neighboring trajectory. This 
deviation from the nominal trajectory is assumed Gaussian with mean 0 so that the sta- 
tistics of the actual trajectory can be completely defined by a six-dimensional Cartesian 
covariance matrix D in position and velocity. For a single off-nominal trajectory, 
however, the actual state of the spacecraft will not be known. It is assumed that the 
spacecraft has been tracked and that a filtering process has been used to estimate or 
predict the actual state of the spacecraft at some time. By the very nature of an esti- 
mate it will be in error. Several factors affect the accuracy of the estimate, such as the 
observability of the spacecraft and the accuracy of the observations. Statistically, the 
error in the estimate can be defined by a six -dimensional Cartesian covariance matrix T 
in position and velocity. In order to define the matrices T and D it is first neces- 
sary to define the various types of trajectory perturbations. The actual state deviation 
AX d and the estimate of this deviation AX e are related by (see fig. 2) 

AX e = AX d + AX t (1) 

where 

AX e = X e - X n 

AX d ^ x a - x n 

AX t s X e - X a 

Here AXt is the error in obtaining the estimate of the trajectory. With these defini- 
tions the matrices D and T can now be defined. As stated earlier D is a six- 
dimensional covariance matrix describing the statistics of the actual deviations from the 
nominal or 

D = E [AX d AX d T ] = E[(x a - X n )(x a - X n ) T _ 
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Figure 2.- Sketch of trajectories. 


Likewise 

or 


T is a six-dimensional matrix describing the statistics of the tracking errors, 


T = E (aX 4 AX t T ) = E 




Since D is the covariance of AX^ and T is the covariance of AX^, a sample 
estimate of the state deviation AX e could be formed by adding a random sample from 
D to a random sample from T (see eq. (1)). However, this implies that AX^ and 
AXt are uncorrelated which is not the case.^ It is shown in reference 2 (and used in 
ref. 3) that for an optimal estimate e(ax 6 AX^) = Jjoj when AX e = AX^ + AX^. That 
is, the error in the estimate is uncorrelated with the estimate, thus the covariance 
between AX^ and AX^- is 


E (AXd AX t T 


) = E ["(■ 


AX P - 



= E 




= -T 


( 2 ) 


1 The assumption that AX^ and AX^ are uncorrelated is commonly used in 
Monte Carlo analysis; however, the assumption is incorrect. To the authors* knowledge 
the correlation between AX^ and AX^ was first derived for an optimal estimate in 
1962 by researchers at Ames Research Center (ref. 3). 
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An expression for the covariance of AX e can now be found. From equations (1) and (2) 


e(aX 6 AX e T ) = E [(AX d + AX t )(AX d + AX t ) T 

= E (AX d AX d T + AX t AX d T + AX d AX t T + AX t AX^ 

= D - T (3) 

( R, 

Therefore, a random sample of the estimate AX e can be generated from D - T. A 
random error in the estimate AXf. can then be generated from T independent of the 

choice of AX e R since E (aX 6 AXj-^) = [oj. The random state deviation vector is 
given by equation (1) as 

AX d R = AX e R - AX t R 

During the Monte Carlo process, the random estimate of the initial orbit is targeted to 
the final orbit which produces a new set of control parameters. These controls are then 

applied to the actual initial orbit (x a R = X n + AX d R ) to establish the final orbit. This 
final orbit, however, will be in error because the control parameters were calculated 
with the estimate and applied to the actual. In other words, the error in the estimate 
introduces errors into the final orbit. 

Also contributing to the errors in the final orbit are execution errors and the uncer- 
tainty in the spacecraft parameters. The control parameters a, fi, 6, 6, t^, v Q 
cannot be applied or executed exactly as computed because the direction of the thrust 
will be slightly different than desired, the thrust vector will rotate faster or slower than 
computed, and so forth. All of these variations will perturb the final orbit. Since the 
spacecraft will be alined to some arbitrary reference system such as a Sun-Canopus sys- 
tem, an attitude maneuver will be performed to establish the computed direction and 
plane of thrust which are defined by a., /3, and 6. As mentioned, this maneuver will 
be in error resulting in improper alinement. Not knowing the actual attitude maneuver, 
it has been assumed that the variation in the direction of thrust is circularly distributed 
about the computed direction and that the error in the plane of thrust is normally dis- 
tributed. Thus, the task of relating these error sources to the variation in a, /3, and 6 
arises. This problem is addressed in the appendix. The errors in 6 and v Q have 
been assumed uncorrelated and normally distributed with mean 0 and variance o%, that is 



A^o ~N(0,oJ o ) 
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The uncertainties in the spacecraft parameters m Q , m, and r, and the gravitational 
constant m also affect the final orbit and must be considered statistically. These 
error sources are also assumed uncorrelated and normally distributed, that is 



Am ~ 


At ~ N^0,a^ 
A/i -N(o,a2) 


The error model for the total burn time t^ has been assumed representative of a closed 
loop system; that is, the thrust engine cutoff is triggered by a velocity counting acceler- 
ometer. The integrating accelerometer senses equivalent velocity pulses and signals the 
engine to cut off at the end of K pulses where each pulse represents v meters per 
second of velocity. There are three sources of error associated with this process. 

First, the calibrated value of v will be in error by a small amount e where 
e ~ N(o,crp). Therefore, the true value of v = v £a j + e. Secondly, the accelerometer 

may be biased a small amount V^, where V^~N^0,ay^j. Finally, the accelerometer 

signals the cutoff on a whole number of counts. Since it goes to the next full count, the 
velocity added will always be greater than desired but not by more than v ca j meters 
per second. The integral of acceleration due to thrust V is defined as 


V = f * a(t)dt = f * 1-r- dt 

J 0 J 0 m Q + mt 



m 0 + mt^ 
, m o / 


(4) 


Thus, the velocity -to -be -gained is a function of the computed control parameters and is 
given by 


V tbg 



+ mt b \ 
m o / 


(5) 


The number of pulses is 


K = 


1000V tbg 

v cal 


+ 1 


where ||^|| denotes the integer part of £. Now, since the accelerometer was assumed 
to be slightly biased and since the calibrated value of v is in error, the thrust engine 
will actually deliver a sensed velocity of 
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In order to model this effect in the program the change in burn time At b corresponding 
to V a is determined from equation (4) as a function of the true values of the engine 
parameters. These true values are calculated by adding random increments from the 
respective distributions to the nominal values of r, m OJ and m. Equation (4) would 
then take the form: 

V ,- L± A T In 

ct • - • 

m + Am 

where At b is the addition in burn time to compensate for Am 0 , Am, and At. Solving 
for At b gives 



Monte Carlo Method 

In general, the Monte Carlo approach to a statistical problem is to sample repeatedly 
the statistical distributions of the independent variables, evaluate the function of these 
variables with the random samples, and accumulate the statistics of the resulting depen- 
dent variables. For complicated functions the Monte Carlo method is often the only means 
by which the statistical nature of the dependent variables can be found. 

The Monte Carlo method is very useful in the error analysis of an orbital transfer 
for the following reasons. The initial orbit is described by two statistical quantities, 

D and T. In addition, the actual control parameters are known only in the statistical 
sense. Hence, the final orbit as a function of these various statistical distributions must 
be described statistically. The problem, then, is to find the statistics of the final orbit 
and also to determine the statistics of the transfer maneuver and the variation in the fuel 
resulting from the errors in the initial orbit. 

The implementation of the Monte Carlo method to the orbital -transfer error analysis 
is outlined in figure 3. The first step is to generate a random estimate of the state devia- 
tions from the covariance matrix [b - t) and a random error in the estimate from the 
covariance matrix T. Next, the actual state vector and the estimate of the state vector 
are formed. During the actual mission, the transfer maneuver will be based on the esti- 
mate of the initial orbit since the actual orbit will not be known. To simulate this condition 
the estimate of the vector state X e is targeted to the final orbit with VITAP (ref. 1). 
This results in a set of up to six control parameters which describe the transfer maneuver. 


(m Q + Am 0 ) + (m + Am)(t b + At b ) 


m o + Am 0 
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Figure 3«- Flow diagram of the Monte Carlo method as applied to the 
orbital transfer error analysis. 






Since these control parameters cannot be executed exactly, they are perturbed by a ran- 
dom sample from their error distributions. The change in the burn time to correct for 
errors in m Q , m, and r is then calculated from equation (7). These perturbed control 
parameters are then applied to the sampled, actual initial orbit to establish the final orbit. 
This entire process is then repeated and another final orbit is established. Finally, after 
many passes through the Monte Carlo loop, the final orbits are processed and the statis- 
tics of the variations in the final orbit are formed. 


The accumulation of the statistics of the final orbit and the associated parameters 
is rather straightforward. If | is a dependent variable and the Monte Carlo process 
has generated N s samples of £ denoted by £i (i = 1, 2, . . ., N s ), then an estimate 
of the mean of £ is 


I = 



and an estimate of the variance is 


or 



(8) 


( 9 ) 


If the dependent variable is an M-dimensional vector, say and the N s samples are 
denoted by (i = 1, 2, . . ., N s ; j = 1, 2, . . ., M) then an estimate of the mean of £3 
is 



and the M by M covariance matrix Cj k is 

i=N s 

^ / \ / — \ 


or 


i=l 


i=N s 




{ ! { i ' {i 


i=l 


( 10 ) 


( 11 ) 
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The mean and variance of a parameter are very useful quantities to describe the 
statistical nature of the parameter. A bar chart (see fig. 4) is another way of represent- 
ing the statistical distribution. Here the range of the parameter is divided into small 
intervals and the number of samples that fall within each interval is counted. The result- 
ing chart describes the statistical distribution and is called a histogram. 



Figure 4.- Bar chart of statistical distribution. 

The question still remains concerning the accuracy of the distribution of | after 
Ns samples. For example, how accurate is the estimate for the mean of |? (See 
eq. (8).) It can be shown that as N s , the number of Monte Carlo samples, increases the 
difference between the estimate and the true mean decreases. Nevertheless, for a given 
number of samples it would be desirable to put an interval around | and be very confi- 
dent that the true mean was within this interval. If the parametric form of the distribu- 
tion of £ were known, the problem would be much easier; however, this is not the case. 
The problem, then, is to define a confidence interval around the estimate for a nonpara- 
metric distribution. These nonparametric confidence intervals of the "central type" are 
well established (refs. 4 and 5). Briefly, they involve solving the cumulative form of the 
binomial density for the upper and lower probability limits of the estimate. This proba- 
bility interval is then mapped into the desired confidence interval about the estimate. 

The Monte Carlo process requires the generation of random vectors from the 6 
by 6 covariance matrices T and jp - t]. This is done in the following way. First 
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the covariance matrix T is diagonalized by transforming it to the principal axis system, 
that is 

diag(\ v \ 2 , . . X 6 ) = $ T T$ 

where the X^ -terms are the eigenvalues of T and the columns of the transformation 

matrix <$ contain the eigenvectors. The standard deviations in the principal axis sys- 

t 

tern are the square roots of the -terms so that a random vector r] is given by 


V R = {r x \\ /2 ,B.2\\ /2 , . . R 6 X^ /2 ) 

where the Rj -terms are independent random numbers from a normal distribution with 
L 

mean 0 and variance 1 . The random vector rj is then rotated back to the original 
coordinate system to form the desired random vector from T, that is 


AXj. R = 

A random sample from the covariance matrix (jD - t) which is a random sample of the 
estimate of state deviations AX e is found in a similar manner. 

The random numbers Rj can be generated in several ways. For this analysis 
they are found by generating two random numbers uj and U2 from a uniform distribu- 
tion defined on the interval [O, l] and evaluating the expression (ref. 6) 

*!=[- 

Frequently in application the covariance matrices T and [b - t) are found to 
be indefinite instead of positive definite. This condition usually results from numerical 
problems in generating the matrices. When the determinant of these matrices is small 
or when the components of the vector are highly correlated, it is likely that numerical 
problems have been encountered in computing these matrices and that they are indefinite. 
If this is the case, then at least one of the eigenvalues is negative and the standard 
deviation ||x^ is meaningless. This problem can be overcome in two ways; the covari- 
ance matrix can be forced to be positive definite or the negative eigenvalue can be set 
equal to 0. The latter was adopted. 


2 In u 1 J 1//2 cos(27m2) 


COMPUTER PROGRAM VEAMCOP 


Program Description 

The computational algorithm for a Monte Carlo error analysis of a minimum -fuel, 
thrusting transfer between two Keplerian orbits has been combined with the computer 
program VITAP (ref. 1) to form the computer program VEAMCOP (Viking Error Analy- 
sis Monte Carlo Program). It contains a main program and seven subroutines in addi- 
tion to VITAP and its 11 subroutines. The entire program is written in FORTRAN IV 
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computer language for the Control Data 6600 computer system and resulted in a field 
length of 73 OOOg. A copy of VEAMCOP can be obtained from COSMIC (Computer Soft- 
ware Management and Information Center), Barrow Hall, University of Georgia, Athens, 
Georgia. 

Various options are available in VEAMCOP. The six control parameters a, 13, 

6, 6 , tj-p v Q may be varied in order to minimize the fuel or they may be fixed at 
specific values. If one or more of these parameters is fixed, then the minimization 
process varies the other controls to find the best transfer trajectory. The transfer tra- 
jectory is also required to satisfy a number of constraints which are outlined in table I 
and fully discussed in reference 1. The option is available to control such parameters 
in the final orbit as the six orbital elements, the radius of periapsis and apoapsis, and 
so forth. In addition, VEAMCOP will perform two modes of targeting. The first mode 
is the normal minimum -fuel transfer outlined previously. The second mode allows for 
the minimum-fuel transfer with the additional constraint that the inclination of the final 
orbit be between an upper and lower bound. 

Several input options for the covariance matrices D and T have been incorpo- 
rated in VEAMCOP and will be discussed with the definition of the program inputs. 

Program Input 

All input to program VEAMCOP is accomplished by means of a FORTRAN name- 
list CASE. Each of the namelist variables is defined in table II. Many of the inputs 
relate to the targeter VITAP and are discussed in reference 1. The remaining param- 
eters are discussed here. 

The statistics of the initial orbit T and D are input through the parameters 
TRACK and DEV, respectively. For convenience three options are available for the input 
of T and D. The desired option is selected by inputing a 0, 1, or 2 in KEY. KEY 
equal 0 implies that the input parameter XNOM contains the Cartesian state of the nomi- 
nal orbit and that the parameters TRACK and DEV contain the full 6 by 6 Cartesian 
covariance matrices at XNOM. All three parameters, TRACK, DEV, AND XNOM are 
expressed in the areocentric coordinate system. An input of 1 in KEY implies that XNOM 
contains the Keplerian orbital elements of the nominal orbit in the areocentric coordinate 
system. For this option TRACK and DEV contain half of the covariance matrices at 
XNOM in the N,V,W-coordinate system. 2 Along the diagonals are input the standard 

^The N,V,W-coordinate system is defined as follows: 

|v| |rxv| [ v | |rxv| 

where R is the Cartesian position vector and V is the Cartesian velocity vector. 
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deviations, and the correlation coefficients are defined in the right off-diagonal elements. 
If this option is exercised, VEAMCOP calculates the full covariance matrices and puts 
variances along the diagonals and covariances in the off-diagonal elements. It then 
rotates the covariance matrices from the N,V,W-coordinate system to the areocentric 
coordinate system. The remaining option, KEY = 2, implies that XNOM contains the 
Keplerian orbital elements of the nominal orbit in the areocentric coordinate system and 
that TRACK and DEV contain half -full matrices at XNOME in an arbitrary Mars -centered, 
inertial, Cartesian coordinate system. The Cartesian state, at which the covariance 
matrices are defined, XNOME, is expressed in the same coordinate system as TRACK 
and DEV. Unlike the second option, this option defines the covariance matrices by input- 
ting variances along the diagonals and covariances in the right off -diagonals. When 

VEAMCOP encounters this option it first fills the left off-diagonal that was not input, 

.<* 

rotates the covariance matrices to the N,V,W-coordinate system using the state defined 
in XNOME , and then transforms the covariance matrices to the areocentric coordinate 
system by using XNOM. These three options are the result of obtaining the covariance 
matrices from different computer programs and are included for the convenience of the 
user. 

Most of the remaining inputs are either self-explanatory or well documented in ref- 
erence 1. The standard deviations of the controls and the spacecraft parameters are input 
through the array CONT as defined in table II. The input parameter NMC defines the num- 
ber of Monte Carlo cases or the number of final orbits which will be generated to form the 
desired statistics. The input quantity MCPRINT is an output option. For each Monte 
Carlo case a one -line summary of the case is automatically output. In addition, the option 
is available to output a more extensive summary for each case. If MCPRINT is greater 
than zero, VEAMCOP will output a block of data for each iteration in the targeting scheme. 
This output is similar to the VITAP printout and has been included in VEAMCOP as a 
diagnostic tool. Another automatic feature of VEAMCOP should be mentioned here. The 
input parameter MAXIT defines the maximum number of iterations allowed for targeting. 

If the targeter reaches this limit before the process has converged to the minimum fuel 
controls, then the iteration process is stopped and that particular Monte Carlo case is not 
considered in the final statistics. However, before proceeding to the next Monte Carlo 
case, VEAMCOP retargets this divergent case from the beginning and automatically out- 
puts a full description of each iteration. If, however, the MCPRINT option was in effect 
during the bad case, then the detail output was accomplished during the targeting cycle 
and there is no need to retarget this case. Finally, the input parameter KOR has been 
included to control the correlation between AX,j and AXj-. An input of 0 for KOR 
implies no correlation whereas an input of 1 implies correlation (see eq. (2)). 

Many of the input options are more readily defined by their function in VEAMCOP. 
For this reason a detailed flow diagram of the main program is presented in figure 5. 
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Figure 5*- Flow diagram of main program. 
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Figure 5.- Continued. 
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Figure 5 • - Continued 






Figure 5«- Concluded. 


Program Output and Illustrative Example 

A reproduction of a computer output from VEAMCOP is presented in table III. The 
illustrative example is an error analysis of the Mars orbital insertion maneuver. The 
spacecraft approaches the planet on a hyperbolic orbit. At an appropriate time, the 
engine is fired to establish an elliptical orbit. The analysis is concerned with establish- 
ing the errors in the ellipse due to variations in the approach hyperbola and in the burn 
maneuver. In addition, the fuel budget for this ensemble of hyperbolas is determined; 
that is, how much fuel should be allocated to the insertion maneuver to insure that a large 
percentage of the off-nominal hyperbolas can be retargeted to the desired ellipse. 

The first output from VEAMCOP is a complete listing of all the parameters input 
through the namelist CASE (table III). For illustrative purposes the covariance matrix 
of tracking uncertainties, TRACK, and the covariance matrix of state deviations, DEV, are 
assumed spherical. The tracking uncertainties in position and velocity are 20 km and 
20 m/sec (lo), respectively, and the state deviations in position and velocity are 200 km 
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and 40 m/sec (lo), respectively. The value of the input parameter KEY indicates that 
these matrices are half -full matrices with standard deviations on the diagonals and cor- 
relation coefficients in the right off -diagonals. However, since the matrices are spheri- 
cal, the off-diagonal elements are 0. KEY also indicates that they are expressed in the 

N, V,W -coordinate system at XNOM which is the array of Keplerian orbital elements of 
the nominal orbit in the areocentric coordinate system. The sixth entry in XNOM reveals 
that the covariance matrices are referenced to the point on the hyperbola that is 60° 
prior to periapsis passage. The input array, CONT, contains the standard deviations of 
the control parameters, the spacecraft parameters, and the velocity counting accelerom- 
eter. The guidance law is defined by the three input arrays: NOPT, KOPT, and AIN. 

The first six entries in NOPT indicate that four control parameters a, 6, t^, and v Q 
are free to vary in the minimization scheme to find the optimum transfer maneuver. 

The remaining two controls, /3 and 6, are held constant at 90° and 0 deg/sec, respec- 
tively (see GS array). Therefore, the direction of thrust remains inertially fixed through- 
out the transfer maneuver at a right ascension of a and a declination of 6 (see fig. 1). 
The next six entries in NOPT indicate that three constraints are imposed on the fi nal 
orbit. The KOPT array and table I show that the three constraints are 1/a, i, and o> 
of the final orbit. These three parameters are constrained to the values in the AIN 
array. Thus, the guidance law consists of four controls to establish a final orbit with 
a = 20 455 km, i = 35°, and oj = 65°. Since the number of constraints is less than the 
number of controls, the fuel required for the maneuver is minimized. Note from NOPT 
that the second, fifth, and sixth constraints were not imposed. KOPT and table I indicate 
that these parameters are the radius of periapsis r p , the longitude of the ascending node 

O, and the true anomaly v in the final orbit at the end of the maneuver. Although these 

target parameters are not constraints, they are nevertheless included as statistical 
parameters of interest. Also, the namelist printout indicates that 100 Monte Carlo cases 
were generated (NMC = 100), the first case was output in full (MCPRINT = 1), and that 
the correlation between and AXt was considered (KOR = 1). 

The next section of output is a full description of the targeting of the nominal state, 
XNOM, to find the nominal controls and the Lagrange multipliers. These nominal values 
are useful as initial guesses in targeting the off-nominal trajectories. Each iteration of 
the minimization scheme is recorded by a selected set of output parameters. The control 
parameters for the first iteration which are the initial guesses input through the GS array 
are given in the first line of printout followed by the initial guesses at the Lagrange multi- 
pliers. These controls correspond to the velocity increment (eq. (4)) labeled "DELTA V." 
Presented next is the orbit which resulted from applying the initial controls to the nomi- 
nal orbit. Since the initial guess at the control parameters was not the optimum set, the 
target parameters are not equal to the desired values. The errors in the target param- 
eters are presented in the next line of output. These errors are then used to calculate 
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the corrections to the control parameters and the Lagrange multipliers which will improve 
the initial guesses. The second iteration follows the same format as the first. The new 
controls are the result of adding the corrections to the previous controls. This process 
continues until the control parameters converge to the optimum set. For the case con- 
sidered here, 34 iterations were necessary to target the nominal trajectory. 

In the data following the history of the minimization scheme the first line gives 
the eigenvalues for the second partial sufficiency check (ref. 1). If these values are all 
greater than 0, then the solution is indeed a local minimum. Negative eigenvalues denote 
that the solution obtained is not a minimum. For the example case, the minimization 
process possesses only one degree of freedom since four controls are varied to satisfy 
three constraints; and, hence, only one eigenvalue exists. Since its value is greater than 
0, we are sure that the nominal solution is a minimum solution. This solution is pre- 
sented on the next line of output followed by the Lagrange multiplier for the nominal case. 
Finally, the TRACK and DEV covariance matrices are output along with their six- 
dimensional eigenvalues. If the correlation between DEV and TRACK is considered, then 
the matrix [DEV-TRACKj is output instead of DEV. 

After the nominal state is targeted, the Monte Carlo process begins. From the 
namelist printout we see that MCPRINT = 1. Therefore, a full description of the first 
Monte Carlo case is output. The first two lines of output data present the random sample 
of the actual state and the random estimate of this state. After seven iterations the 
estimate is retargeted to the desired target parameters. From the second partial check 
it can be seen that the solution is a minimum. This solution is output next. The actual 
initial conic is presented next followed by the actual controls. Notice that the "con- 
trols due to retargeting the estimate" and the actual controls differ. This is due to the 
errors associated with the execution of the computed controls. In other words, the actual 
controls are the result of adding random errors, described by the array CONT, to the 
computed or retargeted controls. The spacecraft parameters are also perturbed to simu- 
late their uncertainty. The actual thrust and gravitational parameters are presented 
next. The "actual controls" are then applied to the "actual initial conic" to obtain the 
"actual target variables" and the "actual final conic." The next two lines of printout are 
the computed or commanded velocity increment (eq. (5)) and the actual or delivered veloc- 
ity increment (eq. (6)). 

The remaining Monte Carlo cases are summarized with a one-line printout. The 
first column is the number of the case and the second column is the number of iterations 
necessary to converge the off-nominal trajectory. The next six columns are the "actual 
controls" and the next six columns contain the "actual final orbit." The last two columns 
are the delivered velocity and the number of positive eigenvalues. 
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The next section of data presents the mean values and the standard deviations of 
the statistical parameters according to equations (8) to (11). The first row is the expected 
or mean values of the computed controls followed by the expected values of the actual 
initial conic. The expected values of the actual final conic and the actual target param- 
eters are presented in the next two rows. The expected values of the commanded and 
delivered velocity increment are next. Following the expected values are the standard 
deviations of the same parameters. As a diagnostic aid the reconstructed TRACK and 
DEV matrices are output. Notice that the mean values of the random vectors generated 
from these covariance matrices are not 0. This is due to the finite number (NMC) of 
samples used to calculate these values. 

The next section of output is devoted to the histograms of the orbital elements of 
the final orbit, the six target parameters named in array KOPT, and the orbital elements 
of the initial orbit. Also included is the histogram of the delivered velocity increment. 
These histograms give a much better indication of the distribution of the random variables 
than do the means and standard deviations alone. The plots can be read in the following 
manner. On the abscissa is plotted the standard deviation of the parameter considered. 

On the ordinate is plotted the number of Monte Carlo cases that fell between the various 
standard -deviation intervals. The strange symbols on the graphs can be interpreted in 
the following manner. If the symbols +--+ appear at the top of a column, the column 
represents the plotted value on the ordinate. If the symbols +..+ appear, the ordinate 
is the plotted value plus a number between 0.1 and 0.49 times the ordinate interval value. 
The symbol ++++ denotes the plotted value plus half of the interval value. The symbol 
+ | |+ corresponds to the plotted value plus a number between 0.51 and 0.99 times the 
ordinate interval value. This interpretation is shown in table IV for interval values of 2, 

4, and 8. 

After the histograms the cumulative density functions are presented. The first is 
the cumulative probability of the delivered velocity increment V a . Here the 100 veloci- 
ties have been ordered and associated with a probability. The second column is an esti- 
mate of the delivered velocity that corresponds to the probability in the first column. 

The 90 -percent confidence interval about this estimate is given in the third and fourth 
columns. For example, the estimate of V a that is equal to or greater than 95 percent 
of the velocities is 1.604 km/sec. This is only an estimate, however, of the true 95th- 
percentile point. Nevertheless, there is 90-percent confidence that the true 95th- 
percentile point is between 1.578 km/sec and 1.667 km/sec. To express this mathemati- 
cally, we define V a * as 

Pr[v a § Vj ^ 0.95 

and state that 

P r jl.578 s V a * ^ 1.667] = 0.90 


27 



where the best estimate of Va is 1.604 km/sec. The final set of data is the cumula- 
tive density functions of the six orbital elements of the final orbit and the six target 
parameters. 


CONCLUDING REMARKS 

A computer program has been developed which performs an error analysis of a 
minimum-fuel, finite -thrust, transfer maneuver between two Keplerian orbits. The 
method of analysis is the Monte Carlo approach where each off-nominal trajectory is 
targeted to a final conic. Basically, the targeting involves the solution of a constrained 
minimization problem by use of constant Lagrange multipliers and the Newton-Raphson 
iteration technique. 

The accuracy of the error analysis, within the modeling assumptions, is limited 
only by the number of Monte Carlo samples generated. In an attempt to model the physi- 
cal problem, the initial and final orbits have been modeled as Keplerian orbits and the 
statistics of the initial orbit have been assumed Gaussian. In addition, the errors in the 
control parameters and physical constants have been modeled as normal random vari- 
ables. Within these modeling assumptions, the analysis can be performed as accurately 
as desired at the expense of computer time. As an indication of the accuracy 90-percent 
confidence intervals are output for some of the statistical quantities. 

The program is applicable to the deboost maneuver (hyperbola to ellipse), orbital 
trim maneuvers (ellipse to ellipse), fly-by maneuvers (hyperbola to hyperbola), escape 
maneuvers (ellipse to hyperbola), and deorbit maneuvers. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 5, 1972. 
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APPENDIX 


ERROR MODEL FOR THRUST ANGLES 

The initial direction of thrust and the plane of thrust for a maneuvering spacecraft 
can be defined by three angles a, |3, and 6. In order to formulate an accurate error 
model for these three angles it is necessary to know the history of the maneuvers which 
alined the spacecraft to the proper attitude. For example, there are four different two- 
axis maneuvers of pitch, roll, or yaw available to establish the proper direction of thrust. 
To establish the proper plane of thrust adds additional alternatives. At first it would 
seem that the maneuver sequence with the smallest sum of rotation angles would minimize 
the total execution error and thus be the desired sequence. However, the spacecraft atti- 
tude control system may be more accurate in one axis than another which would negate 
this choice. More than likely the maneuver sequence will be required to satisfy certain 
mission constraints such as antenna nulls, and so forth. Therefore, it is not possible to 
determine a priori which maneuver sequence will be followed. The only course of action 
is to assume an error model for the angles a, j3, and 5 which approximates the true 
situation. 

In this appendix an approximate error model is formulated for both the constant - 
inertial thrust case and the constant -pitch -rate case. In the constant -inertial thrust case 
it is assumed that the spacecraft is initially alined to an arbitrary reference coordinate 
system and that it is rotated through a two -axis maneuver sequence to aline the thrust 
engine with the line O-F (see fig. 6(a)). It is further assumed that the execution of this 
maneuver is in error, resulting in the thrust engine being alined with the line O-G. As in 
the Gates error model (ref. 7) it is assumed that the engine alinement error is circularly 
distributed about the line O-F. This distribution is modeled by considering i// as a ran- 
dom variable uniformly distributed between 0° and 360° and by taking p as a normal 
random variable, that is p ~N^0,o-p^. The constant -pitch-rate case is similar with the 
addition of a final roll maneuver to aline the pitch axis perpendicular to the plane of 
thrust. This roll maneuver is assumed to be in error an amount y ~ N Thus, 
this appendix is concerned with mapping the above error sources into the thrust angles a, 
/3 , and 6 . 

Consider the geometry of figure 6(a). The point F defines the desired initial direc- 
tion of thrust, and the desired plane of thrust is defined by a and j3. Due to execution 
errors the thrust alinement coincides with point G which is related to point F through the 
random variables if/ and p. The plane of thrust through point G is defined by 2, the 
desired azimuth angle, and the random variable y. It is desired to determine the three 
angles a ', /3 T , and 6' which describe point G. 
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APPENDIX — Continued 


Z 



(a) Points F and G. 

Figure 6.- Geometry of error model. 


A useful coordinate system is defined by P,Q,W (see fig. 6(b)) where P is along 

A A A. 

the radial line to point F, Q is perpendicular to P and in the a,fi plane, and W 
completes the triad. The transformation from the P,Q,W to the X,Y,Z- coordinate sys- 
tem is given by (ref. 8) 


x 


z 



(Al) 


where 


© 1 = 


cos 5 cos a 

- cos j3 sin a sin 6 

cos 6 sin a 

+ cos /3 cos a sin 6 


-sin 6 cos a 

- cos 1 3 sin a cos 6 

-sin 6 sin a 
+ cos /3 cos a cos 6 


sin /3 sin a 


-sin j3 cos a 


sin /3 sin 6 


sin j3 cos 6 


cos j3 
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APPENDIX - Continued 


or 


Q x = -sin 5 cos a - cos j3 sin a cos 6 


Qy = -sin 6 sin a + cos /3 cos a cos 6 ) (A3) 

Q z = sin /3 cos 5 

Now consider the X ,Y ,Z -coordinate system (see fig. 6(b)) where X* is along the radial 
line to point F, Y' is perpendicular to x' and in an easterly direction, and Z' com- 
pletes the triad. Note that Z' is tangent to the meridian making E an azimuth angle. 

^ T T ^ T x ^ ^ 

The transformation from the X ,Y ,Z - to the X,Y,Z -coordinate system can be considered 
a composite of two rotations. The first rotation is about the Z-axis an amount 8, and 
the second rotation is about the new Y-axis an amount -cp. Thus the desired transforma- 
tion is 


cos 8 
sin 8 
0 


-sin 8 
cos 8 
0 


cos cp 

0 

-sin <p 

0 

1 

0 

sin cp 

0 

cos <p 


x' 

y’ 


or 




X' 

y’ 

z' 


(A4) 


where 


and 


cos 6 cos cp 
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sin 8 cos 
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-sin 8 sin 
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cos (p 
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- f x (f x 2 
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) 
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sin 

<P 

= Fz 





0 

f „ 2 

_ 9\l/2 
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= [?x + 

Fy‘ “j 

J 



(A 5) 


If F is alined with the Z-axis, then F x and Fy are 0 and 8 is undefined. In this 
case 8 is set equal to a.. The angle 2 (see fig. 6(b)) is the angle from the vector 
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APPENDIX — Continued 


Z' to the vector Q. The unit vector z' can be expressed in the X,Y,Z-coordinate 
system as 


Zx 


“ 0 " 

Zy 


0 

_Z’ Z _ 


1 


Z' x = -cos 6 sin <p 
Zy = -sin 6 sin </> \ 

Z* = cos <p 

The angle 2 can now be determined by 


(A6) 


cos S — Q • Z — QxZ x + QyZy + Q z Z z (A7a) 

sin S = (q X z') ■ F 

= (QyZ z - Q z Z-)F x + (Q z z x - Q x Z z ^Fy + (Q x Zy - QyZ x )F z (A7b) 

where Q, z' , and F are given by equations (A3), (A6), and (A2), respectively. 

It is now desired to determine the coordinates of point G. The angular displace- 
ment of G relative to point F is given by the angle p, and the azimuth of G is given 
by 1//. Therefore, the coordinates of G in the x' ,Y' ,z’ - coordinate system are 

G x ' = cos p 

Gy* = sin p sin \fs 

G z < = sin p cos xp 

which can be transformed to the X,Y,Z -coordinate system by (eq. (A4)) 


Gx“ 

Gy 


~ G x r 

G y’ 

G Z_ 


- G ^l 


(A8) 


~ Tt TT ~ TT 

Now consider the X ,Y ,Z -coordinate system (see fig. 6(c)) which is defined in the same 
manner as x', y' , Z. The transformation from the X",Y”,Z”- to the X,Y,Z- 
coordinate system is given by equations (A4) and (A5) where the components of F are 
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Z 



(c) Point G. 


Figure 6.- Concluded. 


replaced by the components of G (eq. (A8)). This transformation is defined by 



(A9) 


Now define Q’ (see fig. 6(c)) which is the vector that establishes the new plane of thrust. 
Its components in the x",Y",Z" -coordinate system are given by 


Q\, = o 

x 

q\, = sin(Z + y) 

q’ ,1 = cos(S + y) 
z 


which can be transformed by equation (A9) 
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APPENDIX - Concluded 


Qx 


0 

Q y 

sin(2 + y) 

S z - 


cos(S + y) 


(A10) 


Both the unit vector to the point G (eq. (A8)) and the unit vector in the plane of thrust 
(eq. (A10)) have been expressed in the X,Y,Z -coordinate system. It remains to determine 
a , p' , and 6'. The vector normal to the plane of thrust n is given by n=GxQ' or 

n x = GyQ z ” G z Qy 
n y = G z Q x - G X Q Z 
n z = GxQy - GyQ x 

Since the angle between n and Z is 0* , we have 
cos p' = n • Z = n z 

sin p' - ]jl - cos 2 0' 


where 0° i /3 T = 180°. The node (see fig. 6(c)) is given by 
Z X n 


or 


N = 


Z x n 


n, : 


N x = r 

sin j3 


N v = 


n- 


x 


sin p 


N z = 0 


which yields the angle a' , that is 
sin a ’ = Ny 
cos a 1 = N, 


'x 


If sin p' = 0, the node is undefined. In this case a' + 6' = 8' and it is assumed that 
a' = a and 6' = o' - a' . If this is not the case, the angle 5’ is obtained from 


cos 6 — N • G — NxG x + NyGy 

sin 5 = (N X G) * n = n x NyG z - nyN x G z + n z xG y - NyG x ^ 

Thus, given the three desired thrust angles a, p, and 5 and the execution errors xf/, 
p, y, the resulting thrust angles a' , p\ and 6* can be determined. 
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TABLE I.- TARGET PARAMETER REQUEST KEYS 


Input 

value 

Input parameter 

KOPT (1) 

KOPT (2) 

KOPT (3) 

KOPT (4) 

KOPT (5) 

KOPT (6) 

1 

a 

e 

i 

cu 

0 

V 

2 

1/a 

r a 

a Latitude at 

reference date 

a Longitude at 
reference date 

a Longitude at 
reference date 

a True anomaly at 
reference date 

3 

r a 

r P 

B • T 

a Latitude at 

reference date 

a Latitude at 

reference date 

a Longitude at 
reference date 

■ 

Orbital period 

5 • R 

>x<^ 

k Declination 
of S 

b Right ascension 
of S 


5 

Voo 


x 

c Latitude of 
landing point 

j 

a Sun angle at 
landing point 



a The value of PERJD and REFJD must be input (see table II). 

^S = Incoming hyperbolic asymptote. 
c The value of PER must be input (see table II) . 

^The value of PER, SLAT, and SLON must be input (see table II). 


00 

-a 




TABLE H.- DEFINITION OF INPUT PARAMETERS FOR PROGRAM VEAMCOP 


Program 

symbol 

Mathematical 

symbol 

Dimension 

J Units 

TRACK 

T 

(6,6) 

km and sec 

DEV 

D 

(6,6) 

km and sec 

CONT 

CT p , o y , blank, o^ 

°v 0 ’ °tn 0 ’ °*> a T 

V CT e. °v b 

11 

deg, deg, deg/sec, 
deg, kg, kg/sec, 
kN, kmS/sec^, 
m/sec, m/sec 

XNOM 

X n 

6 

deg, km, and sec 

NOPT 


12 

None 

KOPT 


6 

None 

AIN 


6 

km, sec, and deg 

GS 


6 

deg and sec 

GL 


6 


HP 


6 

deg and sec 

VI 


6 

rad, sec 

NSTEPS 


1 

None 

MASS 

m o 

1 

kg 

DMASS 

m 

1 

kg/sec 

THR 

T 

1 

kN 

PERJD 


1 

days 

REFJD 


1 

days 


Definition 

Rectangular Cartesian covariance matrix of errors 
in the estimate of state (see KEY) 

Rectangular Cartesian covariance matrix of state 
deviations (see KEY) 

Standard deviations associated with control param - 
eters and spacecraft parameter; the last two 
relate to velocity accelerometer (see eq. (6)) 

Initial nominal orbit expressed in either Cartesian 
elements or Keplerian orbital elements (see KEY) 

Integer array denoting free control variables and 
selected target parameters, NOPT (1 to 6) corre- 
sponds to a, /3, 6, 8, t b , v a ; NOPT(K) = 0 
for Kth control fixed and NOPT(K) = 1 for Kth 
control free where K = 1, 2, . . .,6; NOPT 
(7 to 12) corresponds to constraints; NOPT(K+6) * 0 
for Kth constraint omitted; NOPT(K+6) « 1 for Kth 
constraint considered 

Integer array denoting the specific target parameter 
chosen (see table I) 

Array of values for target parameters denoted by KOPT 

Initial values (guesses) of controls or, /3, 6, 8, t b , v 0 
(see fig. 1); these values will vary or remain fixed 
depending on NOPT (1 to 6) 

Initial guesses on Lagrange multipliers; if not input, 

GL(1 to 6) = 1 

Increments of controls for numerical partial derivatives; 
if not input, HP(1 to 6) = 0.6, 0.6, 0.6, 0.006, 1, 0.6 

Maximum allowable step size for controls during Newton - 
Raphson iteration; if not input Vl(l to 6) = 0.30, 0.30, 
0.30, 0.0003, 50, 0.03 

Integer denoting number of segments used for power 
series solution to equation of motion; if not input, 
NSTEPS = 10 

Initial mass of spacecraft 

Mass -flow rate * 

Thrust of spacecraft propulsion system 

Julian date of periapsis passage on initial conic 

Reference Julian date for constraints of longitude, lati- 
tude, and true anomaly at a reference time; PERJD and 
REFJD need not be input if these constraints are not used 
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TABLE n.- DEFINITION OF INPUT PARAMETERS FOR PROGRAM VEAMCOP - Concluded 


Program 

symbol 

Mathematical 

symbol 

Dimension 

Units 

Definition 

PER 


i 

deg 

Angular distance from periapsis to landing point 

SLAT 


i 

deg 

Declination of subsolar point in areocentric equatorial coor- 
dinate system 

SLON 


i 

deg 

Right ascension of subsolar point; input SLAT and SLON 
only if KOPT(5) = 5 

ERR 


i 

None 

Newton-Raphson convergence criteria; if not input, 
ERR = 10' 6 

MAXIT 


i 

None 

Integer denoting maximum number of iterations allowed; if 
not input, MAXIT = 50 

UMARS 

M 

i 

kmVsec^ 

Mars gravitational constant; if not input, UMARS = 42 828.4 

MODE 


i 

None 

Denotes program mode: 

1 — normal forward targeting; 

2 - forward targeting, inclination within bounds; 

3 - not allowed in VEAMCOP 

BOUND 


2 

deg 

Bounds on inclination for MODE = 2; BOUND(l) is lower 
bound 

NMC 

N s 

1 

None 

Integer number of Monte Carlo cases considered 

KEY 


1 

None 

Integer denoting input options for TRACK, DEV, XNOM; 

KEY = 0 implies XNOM is the Cartesian state of nominal 
orbit and DEV and TRACK are full covariance matrices 
at XNOM; XNOM, DEV, and TRACK are in areocentric 
coordinate system; KEY = 1 implies XNOM is Keplerian 
orbital elements of nominal orbit in areocentric and DEV 
and TRACK are half-full matrices (standard deviations on 
diagonal, correlation coefficients in right off diagonal) in 
N,V,W-coordinate system at XNOM; KEY = 2 implies 
XNOM is Keplerian orbital elements of nominal orbit in 
areocentric and DEV and TRACK are half -full matrices 
(variances on diagonal, covariances in right off diagonal) 
in any Cartesian system at XNOME; DEV and TRACK 
may be full matrices 

MCPRINT 


1 

None 

Number of Monte Carlo cases for which iterations are 
output; MCPRINT S 0 

KOR 


1 

None 

Integer denoting correlation; KOR = 0 implies no correla- 
tion between AX<j and AXt, that is e(aX 6 AX e T ) = D + T; 
KOR = 1 implies correlation between AXj and AXt, that is 
E(AXe AXe T ) = D - T; KOR = 2 used only on case following 
KOR = 1 and implies correlation; this option implies that the 
desired covariance matrices have been established in the first 
case and used for the second case; DEV and TRACK are not 
input when KOR = 2 

XNOME 


6 

km and sec 

Cartesian state of XNOM in same coordinate as DEV and TRACK; 
XNOME input only when KEY = 2 

VCAL 

v cal 

1 

m/sec 

Calibrated value of velocity counting accelerometer 
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TABLE HI. - SAMPLE OUTPUT 


*CASE 

TRACK = 0. 2E +C2 » O.C* C.O* 0.0, 0.0, 0*0* 0.0, 0.2E+02, 0.0, 0.0* 

0 . 0 , 0.0, 0.0 * O.G * 0.2E+02, 0.0, 0.0, 0 • 0 » O.G, 0.0, 0.0, 

0* 2E-C1 « 0* 0 1 0.0* C. 0 , 0.0, 0.0, 0.0, 0.2E~01» 0.0, 0.0, 

0.0, 0.0, G.O , G.C , G.2E-01 , 

DEV = 0.2E+C3, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2E+03, 0.0, 0.0, 

0.0, 0.0, 0.0, O.C, 0. 2E +03 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

0. 4E-01 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.4E—01* 0.0, 0.0, 

0.3, 0.0, C.O, C.G, 0.4E-01, 

CONT = 0.5E+CO, 0.0, 0.0, 0.0, 0.2E+00, 0.32069E+02, 0.477E-02, 

0.1E-02, C.14E+01, C.1E-03, 0.28E+00, 

XNOH a -0.426E+04, 0.211E+01, 0.3689E+02, 0.6377E+02, 0.5195E+02, 

-0.6E+02, 


NOPt 

= 

I? f i 1# C, It It It Ot It It Ot Of 


KOPT 


2t 3 f It It It It 


AIN 

= 

0.48E8783249327E-04, O.C, 0.35E+02, 0. 65E + 02 , 

C.O, 0 

GS 

= 

3. 28F+C2 , 0.9E + C2, -0.15£<-02, 0.0, 0.26E+C4, 

-0.7E+02, 

GL 

= 

0.1E+C1, O.C, G.1E+01 , J.lE+01, 0.0, 0.0, 


HP 

= 

0.1E-C1, 0. IE-01, C.1E-01, 0.1E-04, 0.1E+01, 

0.1E-01 

VI 

= 

O.JE + OO, C.3E+CG, 0.3F+03, 0.3E-03, 0.5E + C2. 

0.3E+00 

NSTEPS 

= 

6 , 


MASS 

= 

3.32C695E+C4, 


□ MASS 

= 

-0.4772E+OC, 


THR 

= 

0.13345E+C1, 


PER JO 

= 

I. 


REF J 0 

= 

0.0, 


PER 

= 

1 1 


SLAT 

= 

I, 


SLON 

= 

I. 


ERR 


O 

. 

m 

» 

o 

vT> 


MAXIT 

= 

ICO, 


UMARS 

= 

0.428 284E + 05 , 


MODE 

= 

1 , 


BOUND 

= 

I, I, 


NMC 

= 

100, 


KEY 

= 

1 , 


MCPRINT 

= 

1 , 


KOR 

= 

1. 


XNOME 


I, I, I, I, I, I, 


VCAL 

= 

0.3E-01, 


SEND 
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TABLE IE. - SAMPLE OUTPUT - Continued 


TARGETING NOMINAL STATE TO FIND NOMINAL CCNTPCLS AND MULTIPLIERS 
INITIAL CCMC 

SMA -426C.0C0C ECC 2.113COOC INC 36.890000 PER 63.770000 N00 51.950000 TAN 300.00000 


ITERATION 

CONTROLS 

l 

2 3. roccococoo 

9C.CC3C0CC0C0 

DELTA V 

-15. COCO -3 OC COO 

1. 368064761 

0. 

2600.00000000 

-70.0000000000 

LGR MULT 

UOODCCCOOKC 

C. 

I.OOCCOOOCOOO 

1,00000000000 

0. 

0. 

ORBIT 

2*6 -2. 1 ? T506 ? 

. 8C7302626007 

37.2517375971 

67. 2646287248 

52.57C2261015 

34.3537807809 

ERRORS 

-9.3B9C6ti9Tfct9E-C6 

C . 

2. 251 75755714 

2. 26462872478 

0. 

0. 

CORRECT 

- 3. 6 % 

-4.3 0. 

-50 15 

-4.46E+05 0, 

1.6 .21 

0. 0. 

ITERATION 

CONTPQLS 

? 

H. = 143:086(2 

■ro.oodoccooco 

DELTA V 

-19.2324629253 

1. 334333511 

3. 

255C.COOCOOOO 

-54. 8974861284 

LGR MULT 

2 i ? t*C 2 

C . 

2. 56132875187 

1.21494833450 

0. 

0. 

ORBIT 

575* 2, 1 2 1 3 2 2 A 

.015149 854543 

37.C0645C3720 

69.5468356004 

53.2955414614 

57.2522071066 

ERRORS 

"3, l " ? 2 c7C79C 4£ 5 

C. 

2. QC 6 4 50 27 1 98 

4. 54633560043 

0. 

0. 

CORRECT 

2.9 3. 

-2.2 C. 

-5C -1.4 

-5.2 3E + 04 0. 

1.3 5. 5 1E-C2 

0. 0. 

ITERATION 

CONTROLS 

1 , 

23. EC19-221S' 

n.CCEDCCCOCO 

delta V 

-2C.52515C4628 

1. 301304279 

0. 

2EOO.OOOCOOCO 

-56.2486963898 

LGR MULT 

-43e';-..- = 3sr' 

c . 

2. 852851 91 37 7 

1.27303103563 

0. 

0. 

ORBIT 

B73' , 7.3 '32227 

.915533338663 

36.9255226106 

69. 5111 30384? 

52.2260537405 

54.9287831792 

ERRORS 

-3.15' 3S9B IE •*'--<■ 5 

C . 

2. 52552261C56 

4, 51113036426 

0. 

0. 

CORRECT 

1.7 7 . 

-.’5 0. 

-op -2.8 

-9. 58E+04 0. 

-1.2 6.64E-02 

0. 0. 

I TER AT I CM 
CONTROLS 

22. = 3~ao = = 7?c 

=<•,.- .'to >:o’c9 

DELTA V 

-2C. 6 75 89 69 302 

1,268367594 

3. 

2450. COCCOOOO 

-59.0150151514 

LGR MULT 

-3S-'»C a . Wo 

c , 

2. 68452162872 

1.33644321952 

0. 

0. 

ORBIT 

’75 ’C-o.-RTCcJ 

.516533502714 

26. 8 54 74 5 74C0 

69.5255791404 

51.7872424336 

50.0323209872 

ERRORS 

: :r 5 C -; 3 

C, 

7 . 6 54745 73993 

4.525579I4C4? 

0. 

0. 

CORRECT 

1 » 2) ) • 

.47 0. 

- c 0 -3.8 

-2.52E+05 j. 

-1.5 -3.54E-02 

0. 0. 

IT5RA T ION 
CONTROLS 

; — — — — — — 

* 56 s - *. 1 77C 


DELTA V 

-?7.CC 4 79405?6 

1. 23551431 7 

3. 

2400.CCCCC000 

-62. 8548882749 

LGR MI.ILT 

- 0 t 6 ; 5 ~ . r ' : ? • 

e 

1.:E5143?C1cO 

1. 30130527592 

0, 

0. 

OP 01 T 

•s = i : p . ~ .3 ’ o r- s 

. 515963622117 

?6.ec4<;5 1 :s 7 6 

69.6152513203 

51.7134306248 

41.8697212269 

ERRORS 

-3,19'’5'‘-A’2-o;-.= 

r. 

1* 6C4951 r S760 

4.61525132026 

0. 

0. 

CORRECT 

' , o “• • 

’ . 0 0. 

- C C -fco 7 

-1.15=606 0. 

-I.l -.24 

0. 0. 

ITER AT ION 
CONTROLS 

S 

: 7, - of if : c 1 f 


DELTA V 

-iR*Zr75^1.«l6R 

1. 203335625 

0 . 

2350.00000000 

-69.1322447153 

LGR MULT 


r. 

7 . 49 cc 5S?C £ 'F^cf — 0 2 

1.050 32196737 

0 . 

0. 

OR 3 1 T 

’ 5""2E c 

,83 7-01 = 3700 

2». ”621157304 

69. 8547372703 

51.7651630779 

24.3896349320 

ERR DRS 

- a j 5 “ a c 6 * e-+ 7 78'--C’ : 

C. 

1. 762115 7PC=8 

4. 85473727037. 

0. 

0. 

CORRECT 

“Si 5 0, 

-4,0 

_ r ' 3 

-1.41E+07 0. 

1.6 -2.4 

0. 0. 

ITERATION 

controls 

** ), ' 

“S* ■»: : o 3 : 

delta V 

-?2. C 7?*<+0R257 

U L71522Q95 

2300.CC000000 

-56.0717190758 

LGR mult 

-Id)' =*%/ o : 

r , 

:• * 7 t2 54 £ 5 C 1 - 

-1. 32672928359 

0 . 

0. 

0H3IT 


l.i" 3 05 3 3991 2 

■aA.‘?7c?4^5 94 2 

70. 6683046696 

51.3544846745 

48.3994895493 

ERR IRS 

- ■“ « ) i ' 9 ^ ? E 1 “T * C ' 

C • 

? . 7~ 5941 P 

5.6632O40696E 

0. 

0. 

CORRECT 

;• 2 

.54 C. 

- s r - a . 3 

-3.26E + 06 0. 

c . 8 2.1 

0. 0. 
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TABLE in.- SAMPLE OUTPUT - Continued 


ITFR AT I ' 

cjn t r :l s 

A 

5*. : ■ 3 : doocuo 

delta V 

-2Z«.C?542 1SR91 

!• 14OO O 01«*2 

3. 

225C.OCOCOOOO 

-60. 3695316234 

l OR 

- : 4 2 * 5 p : 

c <« 

I t ,45qcs3i i<*4 

. 779Z7 3739795 

. • 

0. 

')R3I T 

) 3; - tj i.r 

'.,r 5r^-?303^ 

?L,745^91CC74 

70*7201569799 

51.5334826009 

39.9234679204 

ERR IRS 

>6471-:- 

r , 

1 * 7 4-*'2 c • F r 73~ 

5.7?C15>979fi; 

. 4 

0. 

C'lRR' CT 

» ■* * 

- * ' o. 

-b:< -s. a 

-1.39E+-J7 0. 

-1.3 -2.0 

0. 0. 

HERAT in's 

CnN T ?'JL s 

T 

'OODCC 

DELTA V 

-l 3, 2R° 3R OQ t>4i 

1,1DPT63255 

2:o?4ccccoooo 

-66.1856608568 

LOR "I UL T 

- : : : ' ? " 6 * * - r - 1 

c. 

i:.2''c=e37ca’ 

-1,1941=312377 

0. 

0. 

OR 13 1 T 

2 ’ 7 ;■ r - 

,.c : - 3R ^7 324 ) 


70. P4295 c 390 r > 

51.6 7CC 836790 

24.9481617939 

c P R 1 p S 

-3, ;7 • : r ;• . • 

r 0 

’ . 7ir: 

5. 8429595 304-5 

« 

0. 

cnp ??r.- 

■* 4 

S’ 7. 

- O _ ’ 7 

-1. 79=+ )3 C, 

-1.33E+C2 -6.72E + 02 

0. 0. 

IT r° MK N 
C-TNTf r;LS 

1 0 

2 5. 

- * v : •: •" 

delta v 

-3 4^2307391 «cg 

1» 1 l ’77C4E2 

7. 

2707.76371777 

-93.3743947107 

LOR mult 


C 0 

-12?.CRE6'45Cv 

-672. =29594337 

: , 

0, 

ORBIT 

- 7't' ‘ 

’ , r.;7 •: '.-.ai it 

2 6*?SE c 76’ , E7ci 

67.179tE7165? 

5 T .9589686461 

321.278861313 

FRR 1° S 

-A, 1 ”* fi “ • C 1 C “ 3 •" 

r 

1 4 C C PS 7fc ? c 7 C 1 

2. 179657! 0=27 


0. 

COR 3 ?CT 

1 3 > c 

.’2 '. 

-13 ?,7 

-l.l a =+0 5 <J. 

-: c i,i 9E +02 

0. 0. 

ITFR A T IO 

controls 

; j *?: - 

>o 'C?07 ; 

DH TA V 

-14. 11 =3 = c 7 -P°6 

i •» CB2°5501 9 

1 . 

2' 57,7637)777 

-80.6713325747 

LGR MULT 

- :- a T e * 


-1 37. f 7 a. 31 ?°:- 

-553. 597478535 

0 . 

0. 

OR B f T 


’ . ’ = ’ ; 7 ' 9 2 0 7 

3 i. “ 9 = 2“<23.-r c 

46. 9859939977 

f ’ .°C5489391 3 

231.794005995 

£82 IP S 

- J, 7= : -i- : - C 2 

c ■* 

1 . ^H?42Et>P- 


o. 

0. 

COKREC' 

, P 

— i 1 2 

2- 0 

-2.9A c + r8 0. 

fcf 2.00E+02 

9. 0. 

IT E- 1 ATior 
CON T R n LS 

1 ' 

= ‘ . v :gc: 

DELTA V 

-1m, 

1.052475462 

P 

2!07,76371777 

-77.7357223413 

LGR Mill T 


<■, 

-C 7 . t^4-03°6C- 

-352. 6321 85344 

0, 

0, 

Q8-HT 

- ’ 7 ‘ c t a _ z ~ c . u 

.. i 2 :-‘4-:ts5t 

3-4 r l 64 7 7E 3 

66. 788 7 45773 ri 

51.8305616525 

342.726239762 

= PB IF S 

— c f c 4"” , “7"'‘ ' c r : - - 1‘ " 

r , 

1*S1P1C4775P- 

1,73674577697 

■» 

' 4 

0. 

CORR CCT 

• > ' ■ — ' ? • 

-.= 5 7. 

-5C 3c m 

-’,57-+03 C. 

2.<J?f+C2 3.07E+02 

0. 0. 

ITFR A T iC’», 
C. ON T F r L S 

: t 4 r-.: 7 

? r . urcoco 

LElta V 

-3 4. 6 6C2405 75S 

1,0227^456! 

2057,76371777 

-74. 3432 6C629A 

IGF >1IILT 

- ’ ' 3 

f - 

?:.r.7 c 2i'>c7 a " 1 

-47. 1 30147034= 

3. 

0. 

0R3I T 

- 17 : 1 :-3 

.. 1-.71D37'. 2^-S 

3*:. 47335 4* 1 3’ 9 

fc6, 58’l6l336f 

5! .7201630909 

354.423300554 

FRP IRS 

- ; , 1 : - 4 - r - - £ f • - ; f 

c <* 

l»4‘ T ?? r 4:3l° 7 

1. 5831 81 33651 

3 

0. 

COR p FCT 

- t ; 3 - € 

-,06 0. 

- 5C 4,4 

-7 P 1 5r+C3 0. 

7, 37- +C2 4.98E + 02 

0. 0. 

I T L R A T I 
ccintf ns 

’ -r , 1 ^“4'0*<.*: 

3C"i:*o^c 

CELTA V 

-J 5. 1 15 EE4R 150 

.9924952522 

0. 

2007,76371777 

-69.917150381 8 

LGR ‘MIL T 

- l “ ; : ' 1 ' ^ . t : 

c * 

1 *' J7 0 62^35 34i 

44 1.1 90471 RR7 

7, 

0. 

OR B I T 

o'* 

l » L-;i3 3 7 “’3C-' 

7 b # 4 4 4 4f 447 

oo. 35 364R1 3F‘' 1 

51 .5 320284449 

7.95539958658 

FP» IRS 

- ; - : ~ ► - * ' i f ; • - •* r 

' , 

i ^ 42 3 44f>''‘+4Ar> 

1, 35*A4813COP 

?. 

0. 

CORRECT 


- • - <3 

-3C -.1 

-2.65E + 09 C. 

5,CIF+C3 1.17E+03 

0. 0. 

it:? 

CuO T R CL S 

- *»: a---*" 4:"* 

* r . 'C 30 T. ICC “ 

CELT A V 

-1 B.4£6 

, q 6293C765C 

3. 

i 957,76371777 

-60.7976864238 

LGR *1 Ml T 

J ' ,C fi 

c • 

a 1 4 e * 4 l c 4 : 2 :• 3 

1515. 37747C57 


0. 

IP 6 I T 


1 3 ’ -'2 = = 61 3m 

3o. ?E1C5P*. 514 

66. )4j 6421 7} 5 

C C.951C362«75 

29.096 5271838 

FDR ID S 

- 7 , ? 3 “ 1 = - ; =*'="“ 7_-,. =• 

Cn 

? • ?E 1 G581 c 1 36 

1. 34166?1715 a 

•% 

0. 

CUR-" “CT 


c • 

M'T 1. 3 

-3.24C+39 D. 

1.V55+C4 5.72E+02 

0. 0. 


TABLE ID. - SAMPLE OUTPUT - Continued 


ITE-* AT 

CPNTPTLS ‘ 1. ? 

= • .or u-corcc 

PELT A v 

- 19.^692A3P^to 

.99249525 22 

1 . 

2007.76371777 

-59. 3065330289 

LGR MULT Z £ 4 £ m*, ?4 

( > 

1 T‘4C. ->453233 

2187. 25723597 

e 

0 . 

OP 3 T r -4 13 ‘ 3. I ' 

1 . m-24iig- * 

' 6 . T73P6 32 c f * 

06 . 97193263 34 

50 , 5458584494 

35. 473435720? 

err i c . s - ? , ? *» ~ : c c : " •* c : •_ — : ■ 

c« 

1, 2"’ ^ °fc 3 5fc4 

,97103263336- 

0 * 

0 , 

COP-UC’’ . 7 6 

c • 

5C -.77 

-2.93E+09 0. 

I ,518+04 1 .05E + C3 

0 . 0 . 

JT£R AT J CW *7 

COMTRCLS “ !• ss s r ? 5 »" ri ' 7 

9, .„:7I coo? 

C 8 LTA V 

- 1 r • 6 3 1 3 c R c 9 ! Q 

1, 0223245 61 

0 . 

2057.7637,777 

-59. 576608952? 

LGR MULT 

c. 

32771.66 1' 807 

3241.81 5 C 071 4 

r 

- • 

0 . 

ORBIT a-lPe'i 

: . '-"‘24751,44 

36,16472 3 - 7 7-4 

65.3978ogc545 

5C. 2623506264 

37. 3729737867 

£RP ")R S -s’irrc?' -;7 t " ?r -c ~ 

r , 

1 . 1847 2 377742 

. 8 97 Q 0 Q 0 54 54P 


0 . 

CORRECT ,7 s 

“ 5 _ i C • 

5C -.65 

-2.38F+0R 0 # 

1.79E+C4 1 • 57E+P3 

0 . 0 . 

I T E 0 A T 1 1 ! h 13 

CONTROLS B '*« c 3T 1 * - :-=• -7 

P^.vC777:COCJ 

CELT A v 

-20. C419736438 

1.05247548C 

0 . 

’ 2107.76371777 

-60.2279063669 

LGR MULT -^5, 7 74-~97. 7 

( . 

5065 C. c I 90653 

4815. 072656 3C 

0 . 

0 . 

ORBIT -7j?i.4,7.Vc.:>;«- 

1.7c 7". 35435C6 

36.C94ER16973 

0 5, 8251258631 

5C .C14484917 9 

38.5075228011 

ERRORS - = .Bl*:‘’CBARG' c r.-C'' 

r. 

1 , C945 81 64 726 

.82512586310? 

0 . 

0 . 

CORRECT ,5^ 

-.14 0 . 

5C -.oo 

-2.84E+39 0. 

1.97E+04 2.005+03 

0 . 0 . 

ITERATION >7 

CONTROLS B5.-W»?i-3; 

=; .c' 'jpcccco 

DELTA V 

-20. 1778211? 73 

1.C82955C19 

3, 

2157.76371777 

-61.0325132209 

LGR MULT -L7747 bs-p-jp. a 

c ■» 

773C5. 39678'-l 

6817. 98642107 

0 . 

0 . 

ORBIT -t '7'’7e,A<A;ii 

l . G3 R5923C983 

36. 0C427876 7 4 

65.75349? 5542 

4°. 75C07389B8 

39.2946 274360 

ERRORS 7?. =63 ;518215-C- 

c • 

l.CC4?7876T>7 

.753491554214 

7. 

0 . 

CORRECT .62 7. 

-4.2 5E-72 7. 

50 -.a 7 

-2.818+39 0. 

2.C7F+04 2.33E+03 

o« 0 * 

ITERATION 

CONTROLS 35. Ci -'5S:"e:9 

°C.C:39C0C0CC 

-27. 27C3 9 92 575 

U 113770422 

3. 

2207.76371777 

-61. 9011269808 

LGR MULT -’JE'-R;' - "*'-?.? 

c * 

90557, 1595274 

9146. 38367515 

0 . 

0 . 

ORBIT -29 + 87?. • 75C2 1 

1. 016 1? ?53t>2 9 

05.9141740754 

65.6829592045 

49.5832495871 

39, 8924101790 

ERRORS -5, 7279 1 r?i=5’ 

c • 

,914l74r734C7 

.632959204529 

0 . 

0 . 

CORRECT ,55 0. 

-fc*60E-12 0* 

5C -.«9 

-2.7SE+09 0. 

2.11E+04 2.56E + C3 

0 . 0 . 

ITERATION 21 

CONTROLS 36.-945B3955C 

97. OC 70.77O007 

CELTA V 

-20. 3364009613 

1. 144929173 

0 . 

2257.76371777 

-62.7907459550 

LGR MULT -27B:3c?Re3S,2 

C. 

1 1 2C 26. 76] 224 

11701.71 77647 

0 . 

0 . 

ORBIT 637977, C4E7P7 

. 992 A 34 728111 

25.8244067211 

65.6134656991 

49.3900745743 

40.3776294311 

ERRORS -4,73fl414C5?P75-05 

c • 

.324406721148 

, 613465 6 99C 84 

0 . 

0 . 

CORRECT , -+9 C. 

-4,936-72 0. 

50 -.89 

-2.75E+09 0. 

2.11E+04 2.70E + 03 

0 . 0 . 

ITERATION 72 

CONTROLS 36. 9832744943 

90. 30C0 JCOOOD 

CELT A V 

-20.28572C5C31 

1. 176439009 

0 . 

2307.76371777 

-63.6788765732 

LGR MULT -26C7C38ASA3.6 

c. 

1 33C99. 312576 

14397.4993081 

0 . 

0 . 

ORBIT 154741. 03 T 214 

. 969270601834 

25.7350196194 

65.5449248913 

49.2076444201 

40.7921998720 

ERRORS -4.242539265326E-C5 

C. 

.735C1 5615354 

.544924991276 

0 . 

0 . 

CORRECT .43 C. 

-3.88E-02 0. 

50 -.87 

-2.72E+09 0. 

2.08E+04 2.76E + 03 

0 . 0 . 

ITERATION 23 

CONTROLS 37.4123562437 

90.0000000000 

— > DELTA V 

-20.4245647070 

1.208307932 

0 . 

2357.76371777 

-64.5533571893 

LGR MULT -23792637667.5 

C. 

15388C. 058460 

17160.0590592 

0 . 

0 . 

ORBIT 8 7777. 180QG99 

.945648889115 

35.6460186410 

65.4772582 968 

49.0337530736 

41.1606415506 

ERRORS -3. 7495319847 )2E -05 

0 . 

• 646C18641010 

.477258296755 

0 . 

0 . 

CORRECT .38 0. 

-3.24E-02 0. 

50 -.85 

-2.70E+09 0. 

2.03E+04 2. 77E+03 

0 . 0 . 


rf*. 


TABLE III. - SAMPLE OUTPUT - Continued 


I T£ R AT l ONJ 
COMTF HLS 

2 + 

37. - , *7E->6? 8: <? 

Q O.C 7 70009000 

DELTA V 

-? 0. a*6S61’2’9 

1. 240544221 

0. 

2407.76371777 

-65.4075913255 

LOB MULT 

-3l ; > c :' , V>a,'C,p 

c . 

174159. 5494?7 

10927. 2691447 

7, 

0. 

ORBIT 

ei2:-i.s'”7:t 

♦ <9220 !'?27?454 

35.5E743244t>? 

65.410405616° 

48.8666992510 

41.4980220169 

ERRORS 

- * „ 2 5 r i " 2 c 7f . -i. 3 

c. 

. 55740 244522. 3 

.410405616360 

0 . 

0. 

CORRECT 

, ? i : . 

-?,f»6E-)2 0. 

SO -4 S3 

— 2 • 6 tf F + 09 0. 

1.94E+C4 2.72E + C3 

0. 0. 

ITERATION 
cr N T R OL s 

7 3 • ’ ’.'0* '6=5 7 

CD LDC OOC3 3 

DELTA V 

-20.4680806965 

1.273136445 

7. 

2457,76371777 

-66.2380968438 

LGB MULT 

- j TV 7 TO-'-W. 7 

c , 

193787.115855 

2264fc, 7261799 

T. 

0. 

ORBIT 

r = . ^ 

. RST?il ! 51035 

?5.469! 7S5 °47 

65.3443292785 

'-8.7051558760 

41.8140001523 

£RR TRS 

-?.7iu4 7 P~>C6 9T' 5 b-C 7 

r . 

.46917955472 3 

. 344329378489 


0. 

CORRECT 


- 2 , *65-32 c. 

EC -.81 

— 2.6 6t 09 0 • 

1.89E + C4 2. 63E + 03 

0. 0. 

I TSRA Tier; 

controls 

2 6 ----------------- 

a -.crpi: ; 7r o 

77.v;t 'cocci 

DELTA V 

-20.81221 579*2 

2. 306 3 53475 

0. 

25C7. 76371777 

-67.0431633253 

LGR MULT 

-■m; >’a-t • 

c « 

l' 26*2, 1 34732 

25273. 8738856 

0. 

0. 

ORBIT 

H3 C !. 7' ’"-<.'4 

, ?74 ?42?43I72 

35.28133177^3 

65.2790192932 

48.5480 806 51 8 

42.1150567269 

ERRORS 

-:,3e:‘~=; = r .--.3:-r > 

r . 

• ?8! ?€* T7t es7 

. 279019293! 6° 

3. 

0. 

C3P7 ECT 

i \ t 

t -■* t 

-2 ■» 5 3 b - 12 G« 

50 -* 7 8 

-2.64F+C9 0. 

1.80T+C4 2. 5 0E + 03 

0. 0. 

ITERATION 

CONTROLS 

? 

f s-jfCt 

5 ; . :<■ -t \-ooc 7 

CE L r A V 

-?u.f 2P0’25364 

i. 339544501 

7. 

255’. 76371777 

-67.8220668661 

LGP MULT 

•ivisr:”,!, 

f a 

2 5 Cft9, 1 EC4F3 

27770.1366865 

7. 

0. 

ORBIT 


1 c 3 72 87 sps ?1 3 

35.294C90C Q4 

65.2145005565 

4P, 3946581000 

42. 4C583 7 1046 

ERRORS 

» 76'- 53 ST' 7 

C • 

♦ ?9*C »C4FP39’ 

, 214530556B7 C 

0. 

0. 

CORRECT 


0 . 

-7 -.’ = 

-2.6 2E + C'9 0. 

1.71F + C4 2.33E+03 

0. 0. 

ITT RAT IU>, 
CONTROLS 

2 i 

'8. ^c‘7 * ' 

R'.c 7 :oo7o: 

colt a V 

-20. 56395261* 3 

1. 373 •'391745 

3. 

2607.76371777 

-68. 5745503177 

LOR mult 

M 

t * 

24776 4, 668863 

33100, 7463196 

J a 

0. 

cm it 

I 7 5t Ik Z” a ? ~* 

.”26i j385S7<,4 

7 3.20 741 7R17.7 

65,1308572300 

48. 24427C6605 

42.6901077250 

ERRORS 

-I,?* 555 ” V.-7,--( -- 

l • 

. 2 C 7 4 1 781 ’C45 

, l 50853199567 

•3. 

0. 

CORRECT 

' 

-2*048-32 0# 

”0 -.73 

-2 • 6JE +°9 0. 

1.61F + 04 2. 1 3c + 03 

0. 0. 

ITE° AT I CM 
CGN T ROLS 

? 1 . 

; * • : *• ? / a ; ^ : r 

S C • «.'0 ‘O'* " 3 X ? 

DELTA V 

-20.!P0"^P745H 

1. 40 7*46980 

1. 

?A=7. 76^71 777 

-89. srm i^ri qa. 

LG ^ mijlT 

-•t'.B-7RR7R'.7,E 

C « 

26385". 7682=° 

72231 . 1 477C88 

0 . 

0 . 

OR 31 T 

2 ^ : * * , ? * : e i x 

• 8 r 1 ■545 s: 5765J 

25. 121 684224? 

65.0882646869 

48.7565130322 

42.9717922562 

ERRORS 

** 9 C H-C p 

o« 

• I 2! 6P4224J] o 

6. 826^6868694erE-:2 

T. 

0. 

C IP?ECT 

,r: 2 1 

725-7' C. 

r r -* 7 0 

- 2 . 50 c 4 .p 9 p # 

1.5CE+C4 1.895403 

0 . 0 . 

I T -^.ATirN 
CHNTF f lLS 

' - , , ^ 

= '.0 j.-0'7COC7 

delta V 

-2 0 .81781 6 250 

1 . 4421 73544 

'3, 

2707.76371777 

-69.9981035375 

LGP MULT 

-7', p 

c * 

278817*252''* 5 8 

74116.3766238 

. 

0. 

ORTI T 

21 Z~ 2t ZT - t L Z 

.-770:3 22246 1 

7". 037671 1184 

65.027275578 3 

47.o;i34]97 82 

43. 2574541195 

ERROR S 

-E , 5 ' : 4<-5 t o-_ L 3 

c . 

3. 7 a C i 1 1 c 2 t 8COE- r 2 

2.72c557833C2 7T-'12 

7. 

0. 

CORRECT 

* • i ** ' - 1 ? 2 * 

-1 • r. 

Of -.3? 

-l.’ 7 Ei- T9 0 . 

c 3t +C 3 7.37E + 02 

0 . 0 . 


I T E ° A T I ( 1 N 
CON TP 0|. <; 

■ • - c- - 

9*-. joot 

r c LTA V 

-2°.e215481653 

1. 459^3Dfct 6 

). 

2732.44170106 

-70. 3234344897 

LGP MULT 

-■+8Ea ".-S * P. - 

c. 

EB5E-9. 3227?$ 

34E53. 367978? 

0, 

3. 

OR B ’ T 

' c >. * *'7 Tc i 

.7^-4-a a. 991049D 


64. 999601 $C4s 

47.0921K1P837 

43.4221822080 

ERR IRS 

-= , ii ■ ?f 6 3 ; :i j ! ■ 

. c , 

-s.’2C05C52C6?25-74 

9 8495 45 6453 BE- C4 

' • 

9. 

CORRJCT 

. 5’. -: 3 ' . 

f. •’ 13 -35 0. 

- 2 « 71E~ r ? ?^ 7 3F- r ‘ :> 

2.57E+05 0. 

-1.CPE+C2 -24 

0. 0. 




TABLE HI. - SAMPLE OUTPUT - Continued 


I T “ ? AT I ?: : 


— 


Qc^T A V 


l.s-S^l U5fc 



CONTROLS 1 9* ' 7^-o 

“ • 1 ,7 


Z ~ % o'* ,o" or Qr,' 

-2 r -. 631*6*0654 


3 • 

2732.41458892 

-70. 3212006043 

Lf,P M«IL- 

’*’.4 


r e 

24?2 *C. R4*r3 ^ 


3492S. 3557651 

-■ • 

0. 

0R3IT ’<‘4 

! r c - 



■34 t cccC3C ;c;tg 


6* • 9 C 99 £45 f 

47.9826256659 

43.4265164775 

ERR3RS ~--,z:r:-T 

? ~ ; 4 . c - 

- 1 i 

r , 

-4, ?C7 c Cc?5?l *7f> 

3R 

-3.54-.73741026cE- 

C 7 3. 

0. 

CORRECT -~. : V-i T 

'* 


0* 

3.9CE--5 -3.43E- 

37 

-1.57E + 03 0. 

-3.41E-C2 -3.19E-03 

0. 0. 

ITERATION 73 — 




DELTA V 


1.459411676 



CPM T *ULS 

c : t : ? 


c *'» : r " ‘ 0 r 7 C J ? 

-:0.£?1*66C552 


J • 

2732.41461792 

-70. 32120C9472 

LGP MULT — ►?? 

* 


c« 

2P524C.8:9r 4 


3482°. 3625735 

0 • 

0. 

OPRI T *> 9 * 

* ~ c c a 


« 7 :ti J 7 ?R1 &h** 

3 f • ''COC 7 r ‘''C(* ~ 


6 C .L 3O30OOCO! 

47. 3826257868 

43.4265163467 

FRKORS >,***?*-“ 

F r : ~y- 

-1 :> 

L • 

!. R1 P C «S4C? 5**F- 

12 

1.1277734301 986- 

1 r r. 

0. 

CORRECT 



7. 08=-. 1 0. 

-3.C8F-T9 5. ME- 

l L 

-2.5 0* 

1. 957-C5 3.81E-06 

0. 0. 

ITRRATITN 34 


— 


CELIA V 


1.459411 676 



CONTROLS 39. ,70 cr 

r s z r • 


90, OJ 30 0-0.7:0 

-27. < 7146 6075 1 


0. 

277->. 41461792 

-70. 3212C09472 

LGR MULT - ', = 5‘:424 

sj;, i 


C 9 

2^52*0. 303°33 


3+320. 3625773 

0. 

0. 

HRRIT 



• 7£ 4 6 ^7^ R] o4t> 

’t.OC'iC" 


65.CQ3C OCC COC 

47. 9626257863 

43.4265163468 

ER’-nRS a, *7 

~ r -»: 


C * 

3. 


-2.7234341 T5319E- 

12 3. 

0. 

CORRECT -.,'7r-i; 

« 


— ’.Mi -i.l 0. 

-’.515-11 1.57E- 

11 

3. 3 c. 

-1.72 C -C5 -2.38E-06 

0. 0. 

e: in:NVAi u- c p?p sf-rr 

r pa-? 

t I AL 

C H = CK t • G 1 1 7 3 j 8 5 

+ 06 





NOMINAL CTNTS-LC 









AL D H A ?Q.l-T5r} 

IL* 

T£ 

c L.' , r:o;'T DELTA -1'. (71466 

THDOT n. 

TBUPN 2732.4146 

TAN -70.321201 

NOMINAL MULTIPLIES 









- + 35V'4?-4;5.- 

■A 

•. • 



.904918 ?4g69.3 

62 5 7 49 0. 

0. 



TRACK MATRIX 
VC'C.-KCC 
-5. 79‘.64a3t 14? 

7. 25Ai = 2A="-;A 

0 . 

0, 

A, 

dev-track matrix 

3 96 ~C* 9C j f 
-5. riJArCTS"’ 355 ; F-I 
3. c?7'st;5 ; <ai •>*. 

o. 

0. 

n, 

EIGENVALUES CF T R A C K 

A00.L'"'" , 'S-r )■>' AC r nCC-'OCG 430.00C00030C <-.OOC3C0000003E-04 4.03000C03000CE-04 4.000000000000E-04 

EIGENVALUES IF rSV-TPACR 

39600. CCor 3 A‘>CC.:CCT:oC 396CO.C0OCOC-C 1.200000030000E-03 1 .230000 00 OCOCE- 03 1. 200000000000E-03 


-r. A26f ■’7o2379 2 E-U 5. 95044 3E T. 7 3 A- E_— 1 1 0 . 0. 0. 

.oof ;3ro -1. 4744373960625-10 C. 0. 0. 

-:.A)F-'c3?3 c -*AT : -13 39603. 3C0CJC" 7. C. 0. 

C. 0. I .2000 3000000 3 E-0 3 -1. 41 81 5804811 7 E-1 8 5.688479509798E-19 

C. 0. -8.1) 9770°73605E-19 1. 2CCOOOOOOOOCE-C3 -4. 061351.3O5859E-1 8 

r. 0. 3. 247C87433921E-19 -2. 83 94029 17018E-1 8 1. 200000000000E-03 


-1. 1^7AP--A_'C zzt-l 3 

4 • ,.C TCO'O'.O 

-’.:37P'341=-,C3F-12 


3.5 = 43 89 ( 88 + 825-1 3 
-9.59511A332R6PF-13 
400. CCCCCOOCC 


4. OOOCtO 300099E-04 -4. 72719349 37246-19 
-2.706590324535E-19 4, C C OCCOOOCC 00E-C4 

1.082 3 4 2‘t9«,640S-19 -9. 4f 4676390061 E-19 


1. 9961 59 8365 99E- 19 
-1. 353783768620E-18 
4.QOOCOOOOOOOOE-04 


tn 


05 


TABLE HI. - SAMPLE OUTPUT - Continued 


CASE 1 


ACTUAL 


X 4015 

.7041 Y 57 o9 

.3798 Z 365. 

11756 XO -4.1950086 YO -1.C727570 ZD 2.0021508 

ESTIMATE 

X 4035 

.3120 Y 5745 

.1756 Z 36C, 

97606 XD -4.2056749 YD -1.0434784 ZD 2.0037876 

ITERATION 

CONTROLS 

1 

39. ! 789555655 

90.0000000000 

CELT A V 

- 20 . 6 31 46 6C551 

1.459411676 

0 . 

2732.41461792 

-70.3212009472 

LGR MULT 

-485404 24 8 C 5. 1 

0 « 

285240.809916 

34829.3625749 

0 . 

0 . 

ORBIT 

24943.7369857 

,8C555455oll6 

25. 66C88 5020 8 

65.0666501319 

46.9175979162 

46.6059963810 

ERRORS 

-8.757578619772E-06 

C. 

.660885020795 

6.665013192378E-02 

0 . 

0 . 

CORRECT 

.99 0. 

-.14 0. 

50 -1.5 

1.92E+08 0. 

2.18E+04 6.59E+03 

0 . 0 . 

ITERATION 

CONTROLS 

2 

40. C 737453 56 8 

90.0000000000 

DELTA V 

-20.7686045566 

1.494695623 

0 , 

2782.41461792 

-71.7815389208 

LGR MULT 

-48348178646. ! 

c. 

307CC9. 879594 

41424.2823889 

0 . 

0 . 

ORBIT 

22777.50 1859S 

. 79C 597 210245 

25. 3374523055 

65.0394608468 

46.6403191395 

44.791981183 2 

errors 

-4.534831651682E-C6 

C. 

. 337452305463 

3. 946034678228 E— 02 

0 . 

0 . 

CORRECT 

• 67 

-.13 0. 

50 -1.2 

4.32E+07 0. 

1.246+04 4.51E+03 

0 . 0 . 

ITERATION 

CONTROLS 

3 

40.7446108284 

90.0030C0C0C0 

DELTA V 

-20.9024424225 

1*5 30430444 

0 # 

2832.41461792 

-72.9845558329 

LGR MULT 

-483C4937C90. 7 

C. 

319264.455769 

45936.4334420 

0 . 

0 . 

ORBIT 

209C5. 6139178 

.77C 946802157 

35.C452923291 

65,0089130048 

46.3968321153 

43.4036153632 

ERRORS 

-1. 052761171 549E-06 

C. 

4. 52923291286 4E -02 

8. 913004802707E-03 

0. 

0 . 

CORRECT 

3.19E-02 0. 

-3.11E-02 0. 

12 -.19 

-3.46E+07 0. 

4.98E+02 3.61E+02 

0 . 0 . 

ITERATION 

CONTROLS 

4 

40.8263995737 

90.COOOOCCOOO 

-20.5335050436 

1.533779294 

0. 

2844.00450430 

- 73. 1779724210 

LGR MULT 

-43329565191. 1 

C. 

315862.257150 

46 297.6191810 

0. 

0 . 

ORBIT 

20455. 8C53769 

.765137210649 

24.9995217231 

65.0000593599 

46.3453746879 

43.3097770240 

ERRORS 

— 1.934 34 802940 9E-C 9 

0. 

-4.7827634Ce0C2E— 04 

5. 93598897467 0E-05 

0 . 

0 . 

CORRECT 

-1.25F-C3 0. 

6.56E-C5 0. 

-9.98E-03 1.52E-03 

-3.83E+04 0. 

-24 -11 

0 . 0 . 

ITERATION 

CONTROLS 

5 

40.8251511924 

90.0080000000 

DELTA V 

-20.5334434176 

1.5387720 97 

0. 

2843.99452803 

-73.1764505166 

LGR MULT 

- 48 33 96C7 5C6. 7 

C. 

319838.543064 

46285.5875866 

0. 

0 . 

ORBIT 

20455.0012256 

.765759557682 

34.C9S9999917 

64.9999994687 

46.3454869805 

43.3140621955 

errors 

-3. 168219159244E-12 

C. 

-8. 2991 391 5 3 67 9E -09 

-5. 312967914506 E- 07 

0. 

0 . 

CORRECT 

-6,66007 0 . 

-1.96E-C7 0. 

2.21E-05 -2.72E-07 

-2.18E+03 0. 

-1.23E-02 -4.40E-03 

0 . 0 . 

ITERATION 

CONTROLS 

6 

40.8251505268 

9C.00C03CG0C0 

DELTA V 

-20.9334436134 

1.538772113 

0. 

2843.99455017 

-73.1764507883 

LGR MULT 

-49399605686. 3 

C. 

319838.530771 

46286. 5831 875 

0 . 

0 . 

ORBIT 

20454.5559998 

.765759582583 

35. COCCOOOOOO 

65.0000000002 

46.3454870263 

43.3140618135 

ERRORS 

5.152128723651E-16 

C. 

-1.8189894035^68-12 

1.55523594OO32E-I0 

0. 

0 . 

CORRECT 

l. 59F-10 0. 

6.18E-I1 0. 

-3.22E-C9 6.15E-11 

—4. 3 0 • 

3. 30 E -05 6.04E-06 

0 . 0 . 

ITERATION 

CONTROLS 

7 

40. P25150S269 

90.CCJ33000CO 

-20.9334436133 

1. 538772113 

0 . 

2843.99455016 

-73.1764507882 

lgr MULT 

-48339609691.0 

C . 

319e38. 530804 

46286. 5831935 

0 . 

0 . 

ORBIT 

20455. OOCOOCO 

. 765759582586 

35. COOCOCCOOO 

65.0000000000 

46.3454810263 

43.3140618137 

ERRORS 

1.30I0426C6983E-13 

C. 

0 . 

-9.09 49 4701772 9E-13 

0. 

0 . 

CORRECT 

9. 398-12 C. 

-1.59E-11 0. 

-1.85E-11 1.79E-11 

2.0 0 . 

-1.53E-05 -3.45E-06 

0 . 0 . 



TABLE HI 


SAMPLE OUTPUT - Continued 


EIGENVALUES FOR SECOND PARTIAL CHECK 5.38320071E+06 


CONTROLS DUE TO RETARGETING ESTIMATE 

ALPHA AO. 925151 BETA 90.000000 DELTA -20.933444 THDOT 0. 


ACTUAL INITIAL CCNIC 


SMA -4044.8768 

ECC 2.1533582 

INC 

37.571445 

ACTUAL CONTROLS 
ALPHA 40.817443 

BETA SC. 000000 


DELTA -20.891884 

ACTUAL THRUST AND 
HASS 3146.8506 

GRAVITATIONAL PARAMETERS 
DMASS -.47685832 THR 

1.3331260 

ACTUAL TARGET VARIABLES 

4. 992334489928E-05 4791.78453643 

35 

.0111103246 

ACTUAL FINAL CCNIC 
SHA 20030.709 

ECC ,76077809 

INC 

35.011110 

COMMANDED DELTA V 

1.53877211276 



DELIVERED DELTA V 

1.54036650642 




PEP 64.470445 
THDOT 0. 
MU 42829.704 
66.0397161349 
PER 66.039716 


TBURN 2843. 9946 TAN -73.176451 

NOO 51.289223 TAN 300.40966 

TBURN 2795.5240 TAN -73.500759 

46.7689647178 39.4867025182 

NOD 46.768965 TAN 39.486703 



oo 


TABLE HI. - SAMPLE OUTPUT - Continued 


SUMMARY OF REMAINING CASES 


NMC 

NIT 

ALPHA 

BETA DELTA 

THOOT 

TBURN TA 

SMA 

ECC 

INC 

ARGP 

NODE 

TA 

DV 

NPE 

2 

7 

32.56 

90.00 -16.78 

0.00003 

2596.1 -63.77 

202A9. 1 

.77209 

3A.81 

66.28 

52.22 

44.84 

1.378 

1 

3 

7 

37.06 

90.00 -17. A6 

0.00003 

26A5.8 -71.23 

23A63.9 

.81106 

3A. 6A 

64. 42 

50.81 

51.68 

1.393 

1 

4 

a 

50. 31 

90.00 -19.20 

C. 00003 

2808.7 -75.96 

2A572.0 

. 80A5A 

35. A6 

63.39 

42.05 

40.94 

1.500 

1 

5 

6 

37. A6 

90.00 -19.72 

0.00003 

2705.2 -71.69 

18A59. 7 

.75136 

35.33 

65. 58 

51.42 

45.23 

1.449 

1 

6 

8 

AA. 25 

90. 00 -21.87 

0. 30000 

2896. A -73.19 

2C923.3 

.75811 

35. A6 

63.45 

44.2 7 

38.24 

1.576 

1 

■y 

7 

27. 38 

90.00 -19. A9 

C. 00003 

2 807, A -70. A1 

19117. 8 

.7A280 

35.18 

65.59 

55.16 

34.38 

1.534 

1 

8 

7 

2A.31 

90. 30 -15.93 

C. 00000 

27A9. 8 -70.05 

2C771.6 

.76213 

3A.90 

64.92 

57.11 

30.26 

1.449 

1 

9 

7 

27.73 

90.00 -15. 83 

c. 00000 

2710.2 -73.62 

18705. 5 

.75791 

3A. 66 

67.65 

56.04 

36.19 

1.463 

1 

19 

7 

32.91 

90.03 -20.19 

0.00000 

2907.8 -72.07 

18573. 1 

•7A1A3 

3A.96 

66.99 

51.05 

42.85 

1.548 

1 

11 

8 

23.78 

90.00 -16.75 

c. 00000 

270A.1 -7C.96 

20256.7 

. 7592A 

35. A1 

65.56 

58.96 

25.27 

1.443 

1 

12 

7 

30. 39 

90. 00 -16.00 

c. 00000 

2780.6 -69.71 

1996A. 1 

.75987 

35. 2A 

63.38 

52.69 

44,96 

1.474 

l 

13 

8 

23.39 

90.00 -13.90 

r. 00003 

27C8.6 -73.37 

2C325. A 

.77296 

35.18 

64.47 

59.13 

31.01 

1.434 

1 

14 

7 

41.72 

90.00 -21. o2 

C. 90003 

28A3.6 -73.19 

2C958.5 

.76818 

3A. 88 

65.43 

45.15 

41.25 

1.565 

1 

15 

6 

37.12 

90. 00 -17.1 8 

C. 00003 

2680.6 -72.71 

2C231. 7 

.78283 

35.77 

66.33 

48.94 

50.47 

1.422 

1 

16 

9 

53.71 

90. 00 -20.83 

C. 00003 

2892.5 -79.69 

20131. 5 

.76558 

3A.32 

65.33 

43.31 

35.67 

1.587 

1 

17 

6 

32.53 

90.00 -20.37 

c. coooo 

2726. A -69.07 

220A9. 5 

.77738 

3A.67 

63.66 

50.87 

40.40 

1.474 

1 

18 

6 

3A.A2 

90.00 -21. A6 

C. 00003 

2852. 5 -65. 6A 

26A66. 9 

.79678 

35. OA 

61. 80 

47.81 

41.37 

1.521 

1 

13 

7 

31. 53 

90. 00 -19. 37 

0. 00003 

2857. 3 -69.13 

21879.8 

.77082 

3A.85 

64.77 

51.33 

41.75 

1.541 

1 

29 

7 

A3. 22 

90,00 -19. A9 

C. 00003 

2623. A -70.93 

20205. 6 

• 7703A 

3A. 88 

64.47 

48. 15 

44.86 

1.379 

1 

21 

7 

A9.25 

9C.C3 -21.57 

C. 00003 

2838.7 -77.25 

2C9A7.0 

.77279 

35.18 

65.67 

43.39 

37.17 

1.547 

1 

22 

6 

33.93 

90. CO -19.31 

C. 00009 

2675.7 -67.59 

2C867. 1 

. 7 66A9 

3A.55 

65.21 

48.30 

43.61 

1.429 

1 

23 

a 

A6. 91 

50. CO -23. 3A 

C. 00003 

2877.0 -76.32 

22620.1 

.78661 

34.77 

64.09 

43.96 

39.35 

1.589 

1 

2 A 

6 

33.57 

90,00 -18. 66 

C. 90000 

2686.3 -70.31 

25177. 7 

.81703 

3A.99 

62.77 

50.50 

49.98 

1.423 

1 

25 

9 

35.73 

50. CO -16.01 

0.03006 

2518.3 -69. A5 

19A82.1 

• 77A70 

3A.A7 

65.40 

52.47 

49.67 

1.320 

1 

26 

5 

33.65 

90.00 -21.37 

0.00003 

2778. 2 -71. 13 

21C76.5 

.77072 

3A.69 

64.63 

46.77 

43.30 

1.505 

1 

27 

7 

Afl. 5A 

90. CO -20. 1A 

C, 00003 

2732.9 -76.28 

2C63A.7 

.7731 A 

3A.70 

65.51 

AA. 90 

37.76 

1.489 

1 

23 

9 

53.57 

50. 00 -17.83 

C.OOCOO 

2913.0 -82.31 

2CA58.2 

.77992 

3A.93 

64.78 

41.29 

42.51 

1.573 

1 

29 

11 

53.33 

90.00 -13.22 

C.OCOOO 

3083.0 -R0.6A 

21929.7 

.77567 

3A.96 

54.95 

39.42 

37.86 

1.718 

l 

33 

7 

32. 31 

= C.CC -20.01 

C. 00009 

2590.8 -65.98 

19598.2 

•7AA1 A 

3A.81 

66.10 

50.34 

35.96 

1.409 

1 

31 

5 

A A. A A 

50.00 -19.90 

0.00030 

2710.8 -72.19 

21A81. 6 

.77711 

3 A. 75 

65.22 

45. 40 

41. 63 

1.449 

1 

32 

fe 

36. 39 

50. 03 -20.63 

C. 00003 

2789.0 -69.69 

19575.7 

.75202 

3A.A6 

64.76 

48.52 

41.76 

1.522 

1 

33 

7 

28.21 

50.03 -17.59 

C. 00003 

2685.7 -7C.56 

2C6C9. 5 

.76680 

3A. 5A 

64.55 

54. 86 

34.50 

1.458 

l 

3 A 

6 

A2.A1 

50. 00 -23. A3 

C. 00300 

2673. 2 -72.07 

21710.6 

. 7BA00 

3A.A8 

65.09 

46.36 

42.71 

1.413 

l 

35 

6 

AA.2A 

?C. 03 -19.22 

C. 00090 

2813.6 -7A.89 

21A10. 0 

, 78A25 

3A.60 

64.69 

45.71 

46.35 

1.513 

1 

36 

7 

A7. 13 

50, 00 -13.62 

0.00003 

28A6.8 -75.95 

19862.7 

.76936 

35. 1A 

64.49 

45.48 

47.66 

1.507 

1 

3 T 

5 

A?. 33 

50.90 -20.53 

C.OCOOO 

2711.3 -71.66 

2178C.A 

.78287 

3 5. 41 

65.56 

47.59 

43.30 

1.437 

1 

38 

b 

A3. 35 

50.09 -20.53 

C. 90003 

2730.0 -71.73 

1 881 A, 7 

•7A177 

35.33 

64.88 

44.81 

40.32 

1.440 

1 

33 

6 

A3. 37 

50, C 9 — 20.63 

C.C0003 

2666.5 -71. 55 

2C352. A 

.76877 

34.99 

64.91 

45.74 

43. C3 

1.432 

1 

A3 

6 

36. 33 

50. 93 -18. 73 

C.OOCOO 

2735.7 -63.99 

223 A9.7 

.783A2 

35.56 

64.11 

48.78 

47.37 

1.452 

1 

A1 

7 

28.65 

90.09 -itj.83 

0. 00000 

2695.3 -67.77 

15389. 5 

•7A308 

35.66 

64.68 

55.06 

35.25 

1.440 

l 

A2 

7 

23.16 

90. CC -13. Cl 

C.OCOOO 

2702.9 -68.67 

218A1.7 

. 77A00 

3A.50 

65.20 

54.76 

35.35 

1.449 

1 

A3 

g 

27. 63 

90.09 -1A.C6 

C.C09C3 

26AA. 7 -70. 5C 

2C973. 9 

• 7761 A 

35.01 

66.02 

56.44 

36.40 

1.409 

l 

AA 

6 

33. 3A 

50, CO -22.13 

c. coooo 

2737.1 -66.95 

2C2A6.7 

•7A730 

3A. 92 

63.77 

52. 19 

36.41 

1.463 

1 

A5 

8 

Al.Ol 

59.00 — 21.56 

C.OCOOO 

2592.0 -72.13 

1818A.7 

.71925 

35.07 

65.53 

43.76 

41.02 

1.604 

1 

A 6 

7 

26. 37 

90.C9 -18.56 

C. 00059 

2592.2 -65.96 

21557. 3 

•7738A 

35. 42 

64. 68 

49. 94 

43. 75 

1.368 

1 

A' 

6 

35,33 

50. 3-5 - 2 3. 0 A 

C. 90033 

2785. 6 -63.33 

18813. A 

•736A9 

35.32 

65.63 

50.41 

41.44 

1.498 

1 

A 3 

R 

A3. 5’ 

<=0,00 -19.9 7 

C. COOOO 

2919. 1 -75.89 

21578. 3 

.76761 

34.61 

64.69 

42.05 

3 6.67 

1.598 

1 

A3 

6 

39, ’6 

90, 99 -23.23 

C. 90093 

265A.2 -67,63 

21657. A 

.77333 

35.12 

65.98 

48.20 

42.45 

1.403 

1 

53 

6 

A2.A1 

90.06 -2C.51 

C.OOCOO 

2775.7 -69.91 

1 7555. 5 

.72252 

35.29 

65.99 

44.77 

41.36 

1.480 

1 

•"1 

7 

36. S' 5 

°O.CC -13.2a 

C. 00333 

2866.0 -70.61 

1 C A6 5,6 

.75005 

35.11 

66.39 

48.3? 

46.05 

1.554 

1 

- ? 

8 

A3, 37 

39,90 -19. A 7 

C.OCOOO 

2811.2 -77. A5 

IP1A1. 8 

. 7A656 

34.68 

65.11 

4-3.82 

42.46 

1.516 

1 



TABLE HI.- SAMPLE OUTPUT 


Continued 


5 * 

6 

as. n 

90. e: 

-13. 33 

C, C0033 

2627.1 -63.69 

21692.5 

.78233 

34.60 

64.66 

49.48 

45.05 

1.404 

1 

5* 

8 

24.36 

9 C • " 0 

-15. 82 

C .00000 

2651.3 -71.23 

19477. 0 

.75781 

35.73 

64. 20 

57.63 

31.46 

1.397 

1 

5 s ? 

p 

4 7 . 73 

c o » c r 

-16. 37 

c.ocooo 

2776.2 -70.13 

19643. 3 

.75235 

35.62 

65.C9 

45. 44 

45.01 

1.458 

1 

56 

10 

48. 47 

50. ^3 

-19.91 

r. ooooo 

3C38.5 -79.65 

2C282.2 

.76875 

35.62 

64.20 

41.33 

43.83 

1.671 

1 

5? 

10 

45.33 

50. rc 

-23.51 

f. ooooo 

299 7 . 5 -75.66 

1 8020. 9 

.71630 

34.53 

6 5.66 

*0.68 

35.29 

1.667 

1 

58 

8 

33.63 

90. CD 

-18.55 

c.ocooo 

2550.7 -71.27 

19125.0 

.77341 

35.64 

66.49 

49.97 

50.13 

1.317 

1 

5^ 

6 

36. 79 

on.r- 

-16.5 7 

c. ooooo 

2587.1 -70.33 

2C9C3.0 

.78149 

34.98 

66. 13 

49.91 

45. 5£ 

1.402 

1 

6 D 

8 

30. "6 

90. r r 

-16.23 

C. 2CC03 

2545. 6 -69.89 

2C443. 4 

.77363 

35.51 

64.15 

55.38 

37.72 

1.370 

l 

61 

7 

3 3 . 7 6 

on,.*'-. 

-23.44 

C. 00003 

2396.2 -67.30 

24996.0 

.78402 

35.35 

64.33 

44.9 7 

38.02 

1.558 

1 

62 

8 

37.57 

9c. :: 

-20. 6 7 

C. 00007 

2514.3 -63.75 

19633. 6 

.76664 

35.23 

67.56 

50.24 

42.56 

1.326 

1 

68 

e 

42. 3 6 

c 0. 

-?0. 4 7 

C. 00003 

2562.7 -72.62 

22767.4 

.80945 

34.76 

63.55 

50.14 

49.25 

1.338 

1 

6* 

7 

33. 32 

5-%CC 

-23. 98 

C. CC007 

2e61.1 -72.63 

19C50.1 

.74715 

34.97 

66.95 

46.53 

42.54 

1.544 

1 

6 5 

8 

52.43 

50. 0T 

-l 7 . 1 5 

C, 30003 

2616.8 -76.22 

24997. 2 

.82011 

35.47 

63.00 

44.25 

43.57 

1.385 

1 

65 

7 

37.33 

5 C * r C 

-17. 5 7 

C. 00007 

2743.3 -70.40 

19278.0 

.75414 

35.02 

64.38 

55.08 

41.91 

1.456 

1 

67 

7 

26. 36 

90. 0D 

-21.90 

C.OCOOO 

2744.7 -67.47 

18651.6 

.72914 

35.03 

65.37 

52.00 

36.02 

1.442 

1 

68 

5 

33.55 

50. rr 

-21.1 5 

C.COOOO 

2652.9 -71.36 

19931. 2 

.76788 

34.83 

66.57 

49.48 

42.62 

1.417 

1 

63 

s 

33.53 


1 

►-* 

<4J 

. 

o 

0.00000 

2777.0 -70.53 

1E789.6 

.74616 

34.59 

66.36 

51.14 

43.09 

1.482 

1 

7 0 

7 

31.7- 7 

CC.0D 

-!o. P2 

c.ocooo 

2624.9 -68.60 

2C56 7. 1 

.77257 

35.31 

64.56 

53.97 

43.96 

1.364 

1 

71 

5 

41.3 6 

°0. 

-21. C3 

c.ncoi 

2721.9 -70.57 

21603. 6 

.77920 

35.54 

64.09 

46.57 

45.25 

1.425 

1 

7? 

6 

43, 53 

'’O. 0 D 

-16.19 

C.COOOO 

7746.7 -72.93 

19964.9 

.76455 

35.19 

66. 1.6 

45. 14 

44.31 

1.466 

1 

73 

7 

29.13 

C 0 , -- 

-19.65 

C. 0030 3 

2695. 9 -65.08 

20116. 7 

. 74C1 8 

35.46 

64.03 

54.31 

34.33 

1.445 

1 

74 

6 

44. 31 

C C. ^ 0 

-26.73 

C • CC 003 

2832.0 -72. 8 7 

19804. 9 

.74031 

34.81 

64.34 

£5.04 

38.26 

1.520 

1 

75 

6 

36. 33 

50. OC 

-22. tZ 

C. -700 03 

2733.1 -67.97 

22399. 1 

.76929 

35.57 

65.00 

48.40 

35.23 

1.494 

1 

76 

9 

33. 23 

90. 

-16.73 

0.00003 

2462.0 -69.64 

20342.6 

.78793 

35.46 

64.38 

51.35 

50. 74 

1.283 

1 

7** 

9 

23,29 

9C. 00 

-11.77 

r.orooj 

2698. 1 -74.37 

2C614. 3 

.77959 

35.35 

65.47 

60.49 

29.70 

1.435 

1 

73 

a 

32.43 

50. ro 

-20.64 

C.OOCOO 

2609.7 -69,50 

22361. 6 

.79359 

35.44 

63.26 

53.62 

43.67 

1.352 

1 

73 

5 

37.39 

5C. f* 

-13.52 

C. ooooo 

2822. 7 -69.8 7 

19747.7 

.75461 

34.86 

65.43 

46.96 

47.37 

1.488 

1 

93 

6 

43. 77 

5C.00 

-20. C 5 

c. ooooo 

2749.8 -72. 15 

22732.4 

• 7 9789 

35.21 

64.52 

45.83 

42.29 

1.470 

1 

81 

5 

39. 22 

C C. 00 

-21.3 3 

C. 30CC7 

2710.4 -69.74 

19724. 3 

.74700 

34.08 

65.41 

48.28 

36.73 

1.472 

1 

82 

R 

25.56 

9C.9C 

-l 7 . 64 

c.cccoo 

2670.3 -68.14 

19593. 9 

.74553 

35.39 

64.62 

55.99 

30.95 

1.432 

1 

93 

6 

39. 36 

C C -,-T 

-22.66 

C.COOOO 

2322.2 -70.62 

19130.2 

.73386 

34.99 

66.34 

45.76 

36. 07 

1.523 

1 

93 

8 

31. 5 7 

^C. OC 

-20. C3 

c.c-cooo 

2512.3 -65.39 

19182. 2 

.74001 

35.54 

65.02 

50.97 

34.30 

1.363 

1 

8 5 

7 

23, 61 

5C.CC 

-20. 62 

c.ocooo 

2816.7 -66.54 

21238.0 

.75258 

35.25 

64.98 

49.65 

36.51 

1.503 

1 

85 

5 

36.33 

-0.^0 

-20.75 

c.ocooo 

2790.0 -70, CS 

ie589. 7 

.73157 

35.03 

65.87 

45.60 

40.63 

1.494 

1 

87 

c 

43.22 

5C. DO 

-16.21 

c. cocoo 

2695.4 -72.23 

20738.3 

.77311 

35.08 

65.04 

45.65 

44.09 

1.439 

1 

83 

8 

2’. 31 

c c.o: 

-lo.33 

c. 00003 

2756.4 -70.02 

19444.3 

.75394 

33.72 

63.90 

56.61 

39.62 

1.414 

1 

83 

3 

44. 17 

50. 

-1~.23 

c.ocooo 

2584.9 -70. 7 1 

221C9.9 

.79835 

35.15 

63.07 

47.43 

51.92 

1.359 

1 

93 

b 

76.15 

5 C.OC 

-23.33 

c.cccoo 

2933.0 -7C.32 

2C139.7 

.73963 

34.79 

65.22 

45.39 

37.18 

1.603 

1 

9L 

7 

46. 25 

5C.C0 

-16.63 

0.00000 

2783.0 -75.24 

22501.0 

.79004 

34.27 

65.09 

41.69 

4C.99 

1.494 

1 

97 

6 

36.73 

90. O' 

-22.71 

C. C-0000 

2734.5 -69.30 

21 Q 80. 3 

.77469 

34.60 

64.32 

51.15 

37.41 

1.487 

1 

93 

11 

35.25 

50.-- 

-13.94 

C. 00003 

2953. 3 -82.10 

18657. 3 

.74414 

35.11 

65.99 

39.95 

33.59 

1.620 

1 

94 

A 

32.77 

c C.0C 

-21.45 

c. cocoo 

2771.0 -65.53 

22044. 9 

. 76C08 

34.66 

64.04 

49.45 

37.95 

1.481 

1 

95 

3 

37.12 

C C. DC 

-15.63 

C. ooooo 

2528.7 -70.11 

18816. 0 

.75953 

35.26 

65.24 

56.56 

39.36 

1.343 

1 

96 

1C 

36. 14 

c c. c: 

-16.2 5 

c.ocooo 

2470.0 -67,36 

18609. 1 

.75724 

34.29 

65.31 

50.46 

46.45 

1.271 

1 

9 7 

8 

29.19 

9C.*3 

-15.6 0 

c.ocooo 

2902. 2 -67.62 

21668. 5 

.75762 

34.88 

64.95 

52.36 

37.63 

1.579 

1 

99 

6 

3’.C4 

9 C . CD 

-17.52 

C.COOOO 

2700.5 -69.15 

2C627.5 

.76817 

34.49 

64.60 

52.92 

42.83 

1.464 

1 

93 

9 

43,79 

cc, OC 

-2-7.1 7 

C.OOCDO 

2989.4 -76.53 

19C7S.4 

.74001 

34.99 

65.09 

42.15 

39.01 

1.640 

1 

103 

7 

23.43 

CC.O'*' 

-15.75 

c. 003C1 

2 7 93.1 -74.72 

1°916, 3 

.77001 

35.27 

66.42 

57.97 

32.73 

1.486 

1 


£» 

CO 



CJ1 

o 


TABLE IB. - SAMPLE OUTPUT - Continued 


EXPECTEO CONTROLS 

ALPHA 37.5238C7 EETA 9C.OOOOOO 


DELTA -194276359 


THDOT 0. 


TBURN 2745. 

.2364 

TAN 

-71.337526 

EXPECTEO ACTUAL INITIAL CONIC 
SMA -4292.1330 ECC 2.1016092 

INC 

36.712273 

PER 

63.561584 

NOD 

52.278343 

TAN 

299.93776 


EXPECTEO FINAL CCMC 
SMA 2C599.435 ECC .76535022 

INC 

35.011366 

PER 

65.009488 

NOD 

49.039143 

TAN 

40.702225 


EXPECTEO TARGET VARIABLES 

4. 883282027536E-05 4806. 23951373 

35.0113664277 

65 

.0094876144 

49 

.0891428478 

40, 

.7022246230 


EXPECTED COMMANDED DELTA V 1.4697004134 

,3 









EXPECTEO DELIVERED DELTA V 1.46915453468 









STANDARD OEVIATICN FCR CONTROLS 
ALPHA 7.7169245 BETA 0. 


DELTA 2.1436425 

THDOT 0. 


TBURN 120. 

C8024 

TAN 

3.5759588 

STANDARD OEVIATICN FCP ACTUAL INITIAL CONIC 
SMA 234.34169 ECC .10176001 INC 

1. 5678806 

PER 

2. 6831534 

NOD 

2.4226739 

TAN 

2.7828447 


STANDARD OEVIATICN FCP FINAL CCMC 
SMA 1634.7060 ECC 2. C 31 70 389E- 12 

INC 

.35314470 

PER 

1.0133377 

NOD 

4.7674693 

TAN 

5.2362056 


STANDARD OEVIATICN FCR TARGET VARIABLES 
3.649916402086E-06 231.229185452 

• 

393144695122 

1. 

01 3 83 76 841 7 

4. 

76746933388 

5. 

23620558985 



STANDARD OEVIATICN OF COMMANDED VELOCITY 9. 483572707831E-C2 
STANDARD OEVIATICN OF DELIVERED VELOCITY 8. A 1 1 1 6 34173 C 8E-C2 


MEANS OF ERRORS GENERATED FRCM TRACK 
421910919754 -1 . 883 328*38 9? A 


553759243777 


-1.339424052054E-03 -1. C 88SC3800790E-04 -4. 13536981 9293E-04 


RECONSTRUCTED TRACK MATRIX 

5 25. 1 58402145 -1 12. 705625356 

-112.705625356 474. 6 377E ”465 

59.4413249192 27.17 51 3385 81 

3.93 76 78 86 05 8"’E-C2 4. 76 «5 40 9948 8 9E -0 3 

1. 041 71439465 -E-C2 -8. 1 776 4C75 2319 E-02 
-8. 3 9095781 53 77 E-C3 3.C62 7 52 C75395E-0 2 


59.4418249092 
2 7. 11513C8681 
454.100399976 
-6. 4C6S58S57675E-C3 
7, 7 626 5 364 3678 E-C 2 
-2.2139897S2846E-02 


3.93767P860587E-02 
4. 768540994389E-03 
-6.406953967675E-03 
3. 99 8 741 4 94 666 E -04 
-7.253153607464E-06 
-1. 929199 007145E-05 


1. 041714394663 E-02 
8. 17764C7523 1 9E-02 
7. 762653643678E-02 
7. 25 31 53607464 E- 06 
3.5C2916311165E-04 
4. 3 87006 11 2452 E-C 5 


-8. 39095781 5877 E-03 
3.0627520758S5E-02 
-2. 213989792846E— 02 
-1.9291 99007145 E-05 
—4* 38700611 2452E -05 
4. 1 977933298 58E-04 


MEANS OF ERRCRS GENERATED FRCM CEV 
-3C. 247343322’ -22. 53 109C 3603 


-18. 61C 9414120 


-3.081234822774E-03 3. 657384269675E-C3 -9. 39026938 21 67E-04 


RECONSTRUCTED 0= V MATRIX 


39655.2890502 
4543. 6C 90 1262 
-41 52. 86162675 
1 51 9 7458 9407 
-1.493802667=2 
- 1 . 1 464484 = 5 9 “ 


45t3. 60 50 9262 
52179. 2C5C534 
-1730. 6° 6] 71 09 
• 22C 2 92 7? C R06 
7. 223919251 9P6E-02 
-. 778C95754316 


-4152.86162875 
-1730. 69617108 
36461.4665829 
-2.289210011 51CF-C2 
-1.93437725463 
-. 6C850C0577C5 


-. 151974989407 
.220292730906 
-2. 28921001 1510E-02 
1.68B993404792E-03 
1.004327758003E-04 
-8.718423806383E-05 


-1.493 60266752 
7.233919251P86E-02 
-1.93437725463 
l.C 0432 775 8C03E-04 
1.429312247609F-03 
9. 1441 229643 60 E-C5 


-1.14644845557 
-.778095754316 
-.638500057705 
— 8.71842380688 3E-05 
S.144122964360E-05 
1.312832283470E-03 



TABLE III.- SAMPLE OUTPUT - Continued 


* 

* 

20 -* 
* 
* 
* 
♦ 

15 -* 
* 
* 
* 
* 

10 -* 
* 
♦ 
* 
* 

5 -* 
* 
* 
* 
* 

0 -* 


1 

1 

1 

l 

1 

1 

1 

1 

— + 
1 
1 


4—4 
1 1 


1 1 

-4 4 l 4 4 

**_**-**-.**! **i**i**i**i**i **!**]_**i** i**i** i **i **-**-**-**- **-**- 

- 3* * * # *2 * * * -1 ** C* L* 2* * 3 4 


SEMI-KAJOR AXIS 
SECOND CONIC 


TOTAL N 


M E AN 


2C 509* 4353698 


STD lt>34. 70598354 


15 


4 4 

i l 


10 -* 
* 
* 

* 

5 -* 

* 

* 

* 

* 

0 -4 


1 

1 

+ 

L 

l 

\ 


1 

4 4 

I 1 3 1 L I l \ 1 4 — 4 4 

I I II 1111 1111 

**-**_*#;** l #* 1 ** \ #* 1**1 **i**i**i **i**i**i **-**-#*_**-** -.**_**-**. 
* -2 -r -1 r 1 2 3 ^ 


f : C C z N T F ICITY 
SECOND CONIC 


TOTAL N 


N£ AN 


,765^5 32 19 7 5fc 


STD 2.G317C3886C48E-02 


* 100 


= 100 


51 



TABLE III. - SAMPLE OUTPUT - Continued 


15 


10 


4 4 

1 1 


l 

1 

1 

I 

4 4 

1 


* 

0 -* 


1 

1 4 — 4 

4 — 4 1 

1 1 ’ 

1 l 


1 1 


** 1 ** 1 ** 1 *-**, ** 1**1 ** 1 ** 1 ** 1 ** 1 +* 1 **-* A -**-**-**-**- **-**-**_**_ 


-3 


INCLINATION 

Second ccnic 

fcC AN 3.' .01 ns 642 77 STD .3S31AA6S5122 


TOTAL N 



1 1 1 

-4 1 1 


X 

l 

1 

1 

1 

i 

i 

4 4 

1 1 
1 4- 

1 
1 
1 
1 
1 
1 

1 1 
1 1 


1 


»*l**i**l** 1**1**!** 1**2 ** i**i** 1**1 **1**1 **-**-**-**-**-**-**-**- 


-A -3 


-1 C 


l 2 


3 A 


ARGUMENT OF 2ERIAPSIS 
SFCCNC CONIC 

MEAN 65 • OOF AS 76 1 A A STD 1. C 1 3 E 3 76 8 A1 7 


TOTAL N 


= 100 


= 100 


52 



TABLE m.- SAMPLE OUTPUT - Continued 


* 

20 -* 
* 
* 
* 
* 

15 -* 
* 
* 
♦ 
* 

10 -* 
* 


5 -* 
* 


+ — 


1 

1 

1 

1 

♦ — + 
1 1 
3 i 

? 1 

1 


+ — + 

1 1 
1 \ 
l l 
: i 

x +- 
i i 
i 

x 

i 

-+ 
i 
i 

i 

i 
i 

l 

i 

i 


+ — + 
i i 
-+ i 
X 


X 

l 

+ — + 
1 1 

1 1 


0 -* **-**-**1 1 **l**i **1**X#* i**i**i**i**i* *-**_**_**_**_**_**_ **_**.. 

-A — 3 -2 -1 C # 1** # 2 * * * * 3 * # * 


TOTAL N 


LONGITUDE OF ASCENDING NODE 
SECOND CONIC 

^E AN 49*089i42E478 STD 4. 767^693 3 38 8 


25 -* 
* 


20 -* 


+ — + 
l 1 


15 -* 


$ 

* 

10 -* 
* 
* 
* 
* 

5 -* 



TRUE ANOMALY 
SFCGNO CCNIC 


• 7C 2 ??^c 22C 


STD 


5# 2362C5ESS35 


TOTAL N 


* 100 


= 100 


53 



TABLE IE.- SAMPLE OUTPUT - Continued 


20 


15 


* 

— * 
♦ 
* 


5 -* 
* 
* 
* 
* 

0 -* 


+ — + 
1 i 


+• — + 
» 1 


1 


1 

-+ 1 
i 

’ i 
+ — + 


**!#***:< 2**1 **l**i#*l**l**l**l**l**l **-**-**-**-**-*#-**-**-**-**- 

- ” “ - 

♦ -3 -2 -1 C 1 2 3 * 


TOTAL N 


TARGET PARAMETER 
N E AN a. 3fc32S2< 27: 36E-35 


<;m T- f-.ARq 1 £4T 2C R£F— 0* 



TAR w? T PARAMETER 


vp AN 4£C‘3« 2 3952 379 


STO 221.229189452 


TOTAL N 


= 100 


= 100 


54 



TABLE in.- SAMPLE OUTPUT - Continued 


* 

ft 

* 

* 


15 -* 
* 
ft 
ft 
ft 

10 -* 
* 
ft 
ft 
ft 


5 


0 


- ft 
ft 
ft 
ft 
ft 

— ft 


1 1 



**1 **!+*! **i**l**:*M*M**i**i **-**_ **_**_ 




-<• -3 -? -I C 1 


L 


TARGET PARAMETER 3 


P PAN 23.0,1 j.366‘*2'’7 


STD . 3931AA6S5 T22 


TOTAL N 


* 

20 -* 
* 
ft 


15 


* 

10 -* 
♦ 
ft 
ft 
ft 

5 -* 
* 


1111 
\xllll 
+— + 1 1 1 l l 1 

X 1 1 l 1 1 i 1 

1 + — + 12 111 

* 1 1 1 1 1 1 L + — + 

* 1 1 i 1 1 1 i 1 1 ♦— + 

* +— ♦ +— + 1111111111 

- * **!♦* 1* ft 1** 1 *♦1**1 ft ft 1**1 **i**i**i*%i**i** 1 ft ft-* ft- ft ft— ft ft— ft ft— + ft— * ft- ftft- 

-4 -3 -2 -1 C i* 2*** ’ * *4* 


TARGET PARAMETER 4 
MEAN 65. 1C94376144 STD 


!• Cl 38376841 7 


TOTAL N 


* 100 


= 100 


55 



TABLE III.- SAMPLE OUTPUT - Continued 



TARGET PARAMETER 5 


^EAN 4? # -)8E 14264^8 


STD 4*767^6933 38 8 


TOTAL N 


* 

25 -* 
* 
* 
* 

* 

20 -* 
* 
* 
* 
* 

15 -* 


* 

* 

10 -* 
* 
* 
* 
★ 

5 -* 
♦ 


1 

1 

1 

1 

1 

t 

i 

1 

1 

+ — + 
1 1 
1 1 
1 l 
+ l 
1 1 


♦ 

0 -* 


+ — +■ — + 


-+ — + 
1 1 


1 


**_**1**1**’«* ] *«i** 1**1 ** l** 1**1** 1** 1 **-**-**-**-**-**-**-**-**- 
-4 * * -3 * * -2 *1 C 1 2 3 * 


TARGET PARAMETER 6 

MEAN 4C *702 2246 23 0 STD 


5*23620558986 


TOTAL N 


= 100 


= 100 


56 


TABLE in.- SAMPLE OUTPUT - Continued 


♦ 

♦ 

* 

* 

15 -* 
* 
* 
♦ 
* 

10 -* 
* 
* 
* 
* 

5 -* 
♦ 
* 
* 
* 

0 -* 


♦*1**1** 1**1 **i**i**i**i**i*+i**i**i**i* *-**-**-**-**- ♦*-**-* 
- - 

-4 -3 -2 -1 0 1 2 3 A 


SB M [-MAJOR AXIS 
INITIAL CCNIC 

MEAN -4252.13796042 STD 234. 34 16efc47C 


TOTAL N 


* 

* 


■*— 

1 

i 

♦ 

1 





* 


i 

i 

+- 

- + 




* 

+— 

- + 


i 

i 




-* 

1 

i 

i 

i 

i 




* 

1 

i 

l 

i 

i 




* 

X 

i 

i 

i 

i 




* 

1 

i 

i 

i 

1 




* 

l 

i 

i 

i 

i 




-* 

1 

i 

i 

i 

i 




* 

1 

1 

i 

i 

i 




* 

1 

1 

i 

i 

i 




* 

+ — + 

i 

i 

i 

+— 

-♦ 



* 

1 1 

i 

i 

i 

1 

i 



-* 

+— + 1 

i 

i 

i 

i 

i 



* 

l 1 1 

i 

i 

i 

i 

♦- 

- + 


* 

♦ — ♦ 1 i 

i 

i 

i 

i 

i 

i 


* 

1111 

l 

i 

i 

1 

i 

i 


* 

1111 

i 

i 

i 

i 

i 

♦— 

-+ i 

-* 

♦*_**_* *1**1 **1**1 *+1**1** 

i**l **1**1 **1**1 **-**_ **_**_**. 

- 

4 -3 -3 -1 


C 


I 


• . 

2 

3 * 


ECCENTF ICITY 
INITIAL CONIC 

MEAN 2. 101o3917C99 STD • 1C1 76C0C5145 


TOTAL N 


* 100 


* 100 


57 



TABLE in.- SAMPLE OUTPUT - Continued 


* 

* 

♦ 

* 

20 -* 
♦ 
* 
* 
* 

15 -* 
* 

* 

* 

* 

to -* 
* 
* 
* 
* 

5 -* 
* 
* 
* 
* 

3 -* 


1 
1 
1 
1 

— + 
i 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


1 

l 

t 

\ 

i 

i 

i 

i 

i 

i 

«■--+ 

1 ! 

I +- 
1 1 
1 1 


1 

1 

+ — + 

1 ♦ 


♦ c-**-**!**!*^!##!**.^**]**!**!**!**!*#!**!**-^*- *»-**-**-#*-**-**- 


-3 


-1 


INCLINATION 
INITIAL CCNIC 


NF AN 


3c« 71 22732 12 5 


STD 


1 * 58 70 80 5 9 30A 


TOTAL N 


* 

* 

* 

* 

20 -* 
* 
* 
* 
* 

15 -* 
♦ 
* 
* 
* 

10 -* 


* 

♦ 

5 -* 
* 
* 
* 
* 

0 -* 


1 

i 

1 

+ — ♦ 


+ — + 
1 I 


- + 

1 

1 

1 

1 

1 

+ — + 
1 1 


- + — + 

1 ♦ -- + 


-A -3 -2 -1 C I 2 3 A 


ARGUMENT OF PE R I APS I S 
INITIAL CCNIC 

MEAN 63. 561 5838S46 STD 2. 68 31 5 3 35 AS 3 


TOTAL N 


= 100 


= 100 


58 



TABLE III.- SAMPLE OUTPUT - Continued 


* 

* 

* 

* 

20 -* 
* 
* 
* 
* 

15 -* 
♦ 
* 
* 
* 

10 -* 


♦ 

5 -* 
* 


+ — 
1 


+ — + 
1 1 


1 

1 

l 

l 

t 

1 

1 

1 

l 

1 

1 

1 

+ — + 

1 1 


1 
1 
1 
1 
1 
1 

4 — 4 
1 1 


+■ — 4 
1 1 


i 

+ — f 
1 1 


1 

1 

1 

l 

1 

+ — 
1 
1 
1 
l 


1 

1 

+ — 4 
1 1 


0 - * **_**_.**!** 1 ** 1 **!** 1**1 ** 1**1 ** 1**1 **1 **-**-**- **-**_**-**- **-**_ 
-4* " * * - 3* * -2 **-l # C 1 2 3 4 


LONGITUDE CF ASCENDING NCDE 
INITIAL CONIC 

NEAN 52. 278 3425 75 5 STD 2.42267385815 


TOTAL N 


* 

* 

* 


* 




4 - 

- + 



20 -* 




1 

1 



* 




1 

1 



* 




1 

1 



* 




1 

1 



* 



4- 

4 

1 



15 -* 


+ — 

4 

1 

1 



* 

4- 

-4 

X 

1 

1 



* 

1 

i 

1 

1 

1 



* 

i 

i 

1 

1 

4- 

- + 


* 

l 

i 

1 

1 

1 

1 


10 -* 

I 

i 

1 

1 

l 

1 


* 

1 

i 

1 

X 

1 

4- 

-4 

* 

1 

X 

1 

1 

i 

l 

l 

* 

i 

1 

1 

1 

1 

l 

l 

* 

i 

l 

1 

1 

1 

1 

4 4 

5 -* 

i 

1 

1 

1 

l 

1 

l l 

* 

+ — + 

1 

1 

1 

l 

l 

i l 

* 

i i 

1 

1 

1 

l 

l 

l 1 

* 

5 1 

1 

1 

1 

l 

1 

1 l 

* 

+— + 1 

1 

l 

1 

l 

1 

l + — +— + 

0 -* **- 

■ **-** 1**1 **i **i ** 

: 5 : 

A. 

fr*l **1**1 **1**1 **1**1 **-**_ **-**-**. 

-4 

-3 -2 

-1 


c 


i 

2 3 4 


TRUE ANCMALY 
INITIAL CONIC 

FEAN 299. 93^76482? STD 2.78284^72711 


TOTAL N 


= 100 


= 100 


59 


TABLE m.- SAMPLE OUTPUT - Continued 


* 

* 

* 

20 -* 
* 

* 

* 

* 

15 -* 
* 
* 
* 
* 

10 -* 

* 

* 

* 

* 

5 -* 
* 

* 

* 

0 -* 


♦ — + 

1 1 

1 1 

1 1 

1 l 


1 

1 

1 

♦ 4 

1 4—4 

1 1 1 

1 1 4 - 


**_**_ **i**£**)*#i**i**j **i **i **i** i-**x * *1 **-**-**-**-**-**- **_**> 


-2 


-1 


DELTA V 

AN 1.40919456468 


STD 9*46 11624173 C 8F— 02 


TOTAL N 


100 


60 



TABLE m.- SAMPLE OUTPUT - Continued 


CUMULATIVE PROBABILITY 

PROBABILITY 
1.00000003000CE-02 
2.000000COCOOCE-02 
3.0000C0O000OCE-02 
4.00000000000CE-02 
5.00000000000CE-02 
6. 0 000 00000 00 0E—0 2 
7. OOOOOOCOCOOOE— 02 
8. 000000O0C00CE-02 
9. OOOOOOOOCOOCE-C2 
l.OOOOOOOOCOOCE-Cl 
.llOOOOCOCOOG 
.120000000000 
.13000000000C 
. 14000000C ooc 

.i5orooooooor 
.160000000000 
.170000000000 
.180000000000 
.190000000000 
.200000000000 
.210000000000 
. 22Q0CC00C OOC 
• 23OCC0C0C00O 
.24O0C0O0C30C 

• 2 50O0C00CC0C 
.260000000000 
. 2 70GC000C00C 
. 2 800C 0000 000 
.290000000000 
.300000000000 
. 31COGCOOOOOC 
•320OC000C30 0 
. 3 300C0C 10 300 
. 34C0CCC9C90C 
.350000000000 
•36CCCCOOOOQO 
.3700CC03C300 
. 380000000300 
. 390C OOCOOOOO 
.400000 OOOOOC 
.410000000000 
.42000000O03C 
. 4 30000C OC OOC 
.440000000 30C 

• 4 5^00 00 OC 30C 
.4600000 30 000 
. 47CC OCQOC 300 
.48CCGCC30 3GC 
.49300CC0030 r 
.5000 GCCOC 30C 


OF DELIVERED VELOCITY 


VELOCITY 
1.27137620491 
1.28250528350 
1.31722854155 
1.31998490723 
1.32630822293 
1.33830965717 
1.34280545126 
1.352C8127127 
1.35944895499 
1.36261347282 
1.36365835024 
1.36765674797 
1.36955511379 
1.37783224535 
1.37857910265 
1.3B504038137 
1.39276549130 
1.39677603151 
1.40169994270 
1.40333425612 
1.40420347628 
1 .40918756014 
1.40930029980 
1.41330634333 
1.41396 82 0297 
1.41704121815 
1.42164643382 
1 .42292420386 
1.42545279611 
1.42914837215 
1.43226369398 
1.43240492194 
1.43376451763 
1.43450243051 
1.43728686953 
1.43658 74 8586 
1.435 76822459 
1.44036719212 
1.44195186036 
1.44302651200 
1.44499120278 
1. 44894190852 
1.44890624271 
1.-+4S 36134606 
1 .44923762882 
1.45164086535 
1.4560352 4749 
1.45763250311 
1 .45775087315 
1.46251019556 


0 . 

.446641446712 
1.04365184772 
1.27557561650 
1.28239789758 
1.30488595068 
1.31814115267 
1.32025562947 
1.32488931628 
1.33262433186 
1.33962431542 
1.34347474696 
1.35081972788 
1.35699358438 
1.360965975C6 
1.36292809119 
1. 3642C91783C 
1.36772687910 
1. 36923954286 
1.37506371218 
1.37822395235 
1.38109185701 
1. 38704181390 
1.39330596774 
1.38685333756 
1.40121579045 
1.40298511768 
1.4C393225468 
1.40713393194 
1.4C924306411 
1.4109063731? 
1.41351391810 
1.41467605881 
1.41773652414 
1.42174233936 
1.42293344126 
1.42528979156 
1.42867173345 
1. 43167214830 
1.43237003171 
1.43335471402 
1.43424266584 
1.43617513228 
1.43801 108697 
1.43919701695 
1. 440055C791? 
1.441C7215718 
1.44238683722 
1.44378591113 
1.4463456 6492 


1.000000000000+100 
1.33903513505 
1.34810060397 
1.35885593988 
1.36284885340 
1.36575117662 
1.36912334027 
1.37759791523 
1.37973906722 
1.38791217929 
1.39497236692 
1.40030451359 
1.40312 501698 
1.40429012735 
1.40920488651 
1.41042831405 
1.41357219305 
1.41554921954 
1.41989861506 
1.42256645713 
1.42498118430 
1.42378316314 
1.43221204877 
1.43248687652 
1.43386203514 
1.43505702988 
1.43762701935 
1.43896457471 
1.43999153974 
1.44103601687 
1.44245389895 
1.44402512374 
1.44703697893 
1.448B3702000 
1.4989397 6555 
1.44913658712 
1.45081556949 
1.45460866579 
1.45713863088 
1.45771 575732 
1.46113389614 
1.46290925071 
1.463 86877967 
1.46580453615 
1 .469 16C78022 
1.471 82671897 
1.47340218926 
1.47381362744 
1.47768095555 
1.48044287154 



TABLE IB. - SAMPLE OUTPUT - Continued 


.510000C3C300 
.52OOC0O3C00O 
. 5300000 0C00C 
. 543CCC00CC0C 
.550000010 00 o 
.56COOOCOCOOC 
.57COOCOOOOO 0 
. 58CC0003C300 

• 5 90000000 OOP 

• 6COOGOO 3C 00C 
*610000000 OOC 
.620000000000 
.630000000 OOC 
. 6430CCCOC300 
. 6 5 00C COOOQOO 

• 66CC0G00C 000 
.670000000000 
. 680000000000 
. 690000000000 
.700000000 100 
. 710000000000 

• 720000000 OCC 
. 73CCCC00100C 
. 740000009000 
. 750000000000 
. 76 OOCOOOOOOO 
. 77COCOOOOOOO 
• 7800CC000300 
. 79OOCCOOCOO0 
. 80000000000 0 
. 810000000300 

. 820OCCOOO OOC 
. 830000000000 
. 94OOCOCO0OOO 
. 85000000000C 
. B6000C000000 
. 870C0000C30C 
. 88000000C 000 
. 890000000000 
. 900000000000 
. 9100000000CC 
. 92000000000C 
. 93000000000C 
. 940000000000 
. 950000030000 
.960000000000 
. 970C0C03G30C 
.980000000000 
. 99000000000C 
1.00000000000 


1 .46306764924 

1.4641841 3C10 

1. 46644 7466C5 

1.47026344787 

1.47248755124 

1.4733103C535 

1.47381521376 

1.47968186617 

1.480 87498 510 

1.48242400141 

1.48538241715 

1.48653504781 

1.49601585118 

1.48941683627 

1.45389972334 

1.49398853721 

1.49449124487 

1.4983C627C80 

1.45981586052 

1.5033244 8310 

1.50528312196 

1.50687479144 

1.51342507061 

1.515990C6383 

1. 52C30057312 

1 • 520 8C 52677 5 

1.521 9 5 3405 C4 

1.523O054C574 

1.53402427423 

1 . 5403 565C642 

1.54094328666 

1.54391745256 

1.54689800856 

1.54828297771 

1.55387378353 

1,55768526719 

1.56530526584 

1.57286612034 

1.57560409735 

1.57876132093 

1.58726195110 

1.58858544401 

1.59781554028 

1.60329983805 

1.60427455565 

1.61967171707 

1.63954621657 

1.66677445666 

1 • 6 7C 74487504 

1.71 630680926 


1.44879795718 
1.44892405129 
1. 44904658899 
1.44995168567 
1.45291C69669 
1.45648915387 
1.45766622379 
1.45910321810 
1.46267139984 
1.46339939455 
1.46488438999 
1.46768633503 
1. 47102722758 
1.47257130681 
1.47381222767 
1.476 28 64216C 
1.48022422756 
1. 4816 367C 849 
1.48426676571 
1.48625955865 
1.48746375811 
1.4885654C8C3 
1.49272731517 
1.49357084743 
1.49442481047 
1. 498C7630623 
1.49987424769 
1.50351304225 
1.50557556680 
1. 5087 C 872 887 
1.51439856495 
1. 518083C8006 
1.5 20 6C 270730 
1.52163011469 
1.52264711455 
1.53383277460 
1.54044035089 
1. 54176616097 
1.5452579C146 
1.54776811753 
1.55287009122 
1. 5579C099351 
1.56722745290 
1.57423926465 
1.57805047740 
1.58736847729 
1.59254648922 
1.60241484718 
1.60960259948 
0 . 


1.48182530070 
1.48445340429 
1.48624471460 
1.48730667223 
1.48869360772 
1.49140353371 
1.49393518387 
1.49416550348 
1.49564115393 
1.49867974248 
1.50048227000 
1.50357709684 
1.50538572821 
1.50686897679 
1.51293088754 
1.51560279572 
1.51930797355 
1.52064710908 
1.52149310233 
1.52248951573 
1.52754794185 
1.536C0323244 
1. 540 4861 *74 82 
1.54123572681 
1.54387075817 
1.54649633697 
1.54792415792 
1.55169955121 
1.55568624486 
1.56022939964 
1. 56671013745 
1.57295081408 
1. 57524487589 
1.57781101291 
1.58318712240 
1 . 58770302128 
1.58984076881 
1.59714954377 
1.60169612197 
1.60376158649 
1.60772957590 
1.61901251365 
1.63296599211 
1.64921689543 
1.66681280736 
1.66924670492 
1,67920987569 
1.70140202717 
1.71595968099 
1.000000000000+100 



TABLE ID.- SAMPLE OUTPUT - Continued 


PRCR 

SMA 

ECO 

I8C 

ARGP 

NODE 

TA 

1 

2 

3 

4 

5 

6 

• 01C0 

17959.53 

.7163 

33.7206 

61. 8031 

39.4157 

25.2653 

3.778E-05 

4.314E+03 

3.372E+01 

6.1 80E+01 

3. 942E+01 

2.527E+01 

.0200 

18020.91 

.7193 

34.0783 

62.7742 

39.9463 

29.7020 

3.972E-05 

4. 333E+C3 

3. 408E+01 

6.277E+01 

3.995E+01 

2.970E+01 

.0300 

18141.82 

.7225 

34.2657 

62. 9°96 

40.6845 

30.2555 

4. 000E-05 

4. 333E+03 

3.427E+01 

6. 300E+01 

4.068E+01 

3.026E+01 

.0400 

18164.68 

.7291 

34.2851 

63. 0726 

41.2947 

30.9538 

4. 001E-05 

4.389E+C3 

3. 429E+01 

6.307E+01 

4.129E+01 

3.095E+01 

.0500 

18459.65 

.7316 

34.3170 

63. 2624 

41.3274 

31.0130 

4.070E-05 

4.405E+03 

3.432E+01 

6. 326E+01 

4.133E+01 

3. 101E+01 

.0600 

i 9573, 1C 

. 7339 

34.4621 

63.3762 

41.6851 

31.4585 

4.262E-05 

4.433E+C3 

3. 446E+01 

6.338E+01 

4.169E+01 

3.146E+01 

.0700 

18569.65 

.7365 

34.4729 

63.3917 

42.0453 

32.7321 

4. 392E-05 

4.458E+03 

3. 447E+01 

6.339E+01 

4.205E+01 

3.273E+01 

.0800 

I860 9.15 

.7396 

34.4/67 

63.4466 

42.0468 

33. 5890 

4.399E-05 

4.4S7E+03 

3.448E+01 

6.345E+01 

4.205E+01 

3.359E+01 

.0900 

1 96 51. 58 

.7400 

34.4986 

63.5518 

42. 1527 

34.2980 

4.421E-05 

4.503E+03 

3.449F+01 

6.355E+01 

4.215E+01 

3.430E+01 

.1000 

19657.27 

.7400 

34. 5007 

63. 6644 

43.3115 

34.3814 

4. 444E-05 

4.518E+03 

3.450E+01 

6.366E+01 

4.331E+01 

3.438E+01 

.1100 

137C5.50 

.7402 

34.5326 

63.7747 

43.3897 

34.3323 

4.464E-05 

4.525E+03 

3.453E + 01 

6. 377E+01 

4. 339E+01 

3.438E+01 

.1200 

18769.64 

.7414 

34. 5375 

63.9016 

43. 7581 

34. 5038 

4. 472E-05 

4.528E+C3 

3.454E+01 

6.390E+01 

4.376E+01 

3.450E+01 

.1300 

19813.35 

.7413 

34.5549 

64.0293 

43.8224 

35.2253 

4.515E-05 

4.544E+C3 

3.455E+01 

6. 403E+01 

4.382E+01 

3.523E+01 

.1400 

18814.71 

.7428 

34.5943 

64.C4C6 

43. 96C4 

35.2507 

4.523E-05 

4.567E+03 

3.459E+01 

6.404E+01 

4.396E+01 

3.525E+01 

.1500 

1991 5.98 

.7431 

34.5965 

64.0910 

44.2549 

35. 2886 

4. 535E-05 

4. 581E + C3 

3.460E+01 

6. 409E+01 

4.425E+01 

3.5 29E+01 

.1600 

19050.14 

.7441 

34. 5994 

64. C931 

44.2747 

35.3494 

4.536E-05 

4.581E+C3 

3.460E+01 

6.409E+01 

4.427E+01 

3.535E+01 

. 17G0 

19079.30 

.7441 

34.6301 

64.1059 

44.7738 

35.6721 

4.550E-05 

4. 582E+C3 

3. 460E+01 

6.411E+01 

4.477E+01 

3.567E+01 

.1800 

19117.82 

.7455 

34. 6353 

64.1466 

44.8114 

35.9644 

4.560E-05 

4.590E+C3 

3.461E+01 

6.415E+01 

4.481E+01 

3.596E+01 

. 1903 

19125.04 

.7462 

34.6388 

64.2C13 

44.9004 

36.0192 

4. 570E-05 

4. 598E+C3 

3. 464E+01 

6.420E+01 

4.490E+0 1 

3.602E+01 

. 20C0 

19133.22 

. 7466 

34. 6573 

64.2041 

44. 9747 

36.0746 

4.578E-05 

4.6C7E+03 

3.466E+01 

6.420E+01 

4.457E+01 

3. 607E+01 

.2100 

13162.22 

.74 70 

34.6627 

64.3187 

45.0437 

36. 18 80 

4. 591E-05 

4.615E+C3 

3, 466E+01 

6.432E+01 

4.504E+01 

3.619E+01 

.2200 

13277.96 

.7471 

34.6578 

64.3326 

45.1450 

36. 4011 

4. 606E-05 

4.615E+C3 

3.467E+01 

6.433E+01 

4.514E+C1 

3. 640E+01 

.2300 

19369.40 

.7473 

34.6822 

64. 34C 3 

45. 1467 

36. 4C59 

4, 609E-05 

4.616E+C3 

3.463E+01 

6.434E+01 

4.515E+01 

3.641E+01 

.2400 

19444.33 

.7493 

34.6904 

64.3773 

45. 3878 

36.5147 

4.610E-05 

4.61 5E+03 

3.469E + 01 

6.438E+01 

4. 539E+01 

3.651E+01 

.2530 

19465.63 

.7500 

34.6981 

64.3757 

45.4015 

36. 6687 

4. 615E-05 

4.626E+C3 

3.470E+01 

6.438E+01 

4. 540E+01 

3.667E+01 

.2600 

13476.00 

.7514 

04.751? 

64.4163 

45.4401 

36.7262 

4.629E-05 

4.623E+03 

3. 475E+01 

6. 442E+01 

4. 544E+01 

3.673E+01 

.2700 

19462.10 

. 7520 

34. 7554 

64. 4669 

45,4814 

37.1729 

4.634E-05 

4.64CE+03 

3.476E+01 

6.447E+01 

4.548E+01 

3.717E+01 

.2800 

19593.92 

.7524 

34.7686 

64.4703 

45, 6044 

37. 1763 

4. 639E-05 

4.677E+C3 

3.477E + 01 

6. 447E+01 

4, 560E+01 

3.718E+01 

.2900 

19558.24 

.9526 

■>4. 7985 

04. 4879 

45.6471 

37.4084 

4.655E-05 

4.681 E+C3 

3.479E + 01 

6.449E+01 

4.565E+01 

3.741E+01 

.3003 

19633.63 

. 7539 

34.3353 

64.5250 

45.7117 

37.6298 

4.671E-05 

4.689E + C3 

3.481E + 01 

6.452E + 01 

4. 571E+01 

3.763E+01 

.3100 

19642.34 

. 7541 

34. 3383 

64.5532 

45.7422 

37. 7200 

4.709E-05 

4.690E+C3 

3.481E+01 

6. 455E+01 

4. 574E+01 

3. 772E+01 

.3200 

13734.27 

.7546 

34.8393 

64.5622 

45.7620 

37.7647 

4.745E-05 

4. 6 55E + C3 

3.481E + 01 

6.456E+01 

4.576E+01 

3.776E+01 

.3300 

19747. 73 

.7572 

34. 8307 

64.5977 

45.8350 

37.8615 

4.763E-05 

4. 701E+C3 

3.483E+01 

6.460E+01 

4.583E+G1 

3.786E+01 

.3400 

193C4.92 

.9=76 

34.8479 

64.6169 

46. 3642 

37,9523 

4. 771 E- 05 

4.7C5E+03 

3. 485E+01 

6.462E+01 

4.636E+C1 

3.795E+01 

. 3500 

19362. 73 

.7579 

34. e625 

64.6267 

46.5345 

38.0173 

4.774E-05 

4.715E+03 

3.486E+01 

6.463E+01 

4.653E+01 

3. 802E+01 

.3600 

1991 6.35 

.7579 

34. 8762 

64.6618 

46. 57C7 

38. 24C8 

4, 779E-05 

4.717E+03 

3. 488E+01 

6.466E+01 

4.657E+01 

3.824E+01 

. 3700 

19921.22 

. 7581 

34. 8918 

64.6829 

46.7690 

38. 2553 

4.784E-05 

4.719E+03 

3.488E+01 

6. 463E+01 

4.677E+01 

3.826E+01 

.3800 

19964. 14 

.7592 

34.8834 

64.6840 

46. 7725 

39. 0117 

4. 792E-05 

4.722F+C3 

3.488E+01 

6.468E+01 

4.677E+C1 

3.901E+01 

. 3900 

19964. 93 

. 7 5 95 

34.9302 

64. 6850 

46.9569 

39.3506 

4.814E-05 

4.724E+03 

3. 490E+01 

6. 469E+01 

4. 696E+01 

3.935E+01 

• 40CC 

19975.05 

. 7599 

34.9215 

64.6941 

47.4345 

39. 3551 

4. 822E-05 

4.729E+03 

3.492 E + Ol 

6.469E+01 

4.743E+01 

3.936E+01 

.4100 

33020. 71 

.7601 

34.9317 

64.7579 

47. 5944 

39.4867 

4. 846E-05 

4.74CE+03 

3. 493E + 0 1 

6.476E + 01 

4. 759E+01 

3.949E+01 

.4200 

3311 6.71 

.76 38 

34.9608 

64. 7652 

47.8111 

39,6218 

4. 848E-05 

4.759E+C3 

3.4C6E + 01 

6.477E+01 

4.781E+C1 

3.962E+01 

.4300 

23131.47 

.76 21 

34. 9644 

64.7785 

48.1456 

40.3164 

4. 851E-05 

4.77CE+C3 

3. 496E+01 

6, 478E+01 

4.815E+01 

4.032E+01 

.4400 

23139.67 

• 76t5 

34.9724 

64. e7S4 

48.1964 

40.3 994 

4.852E-05 

4.770E+03 

3.497E+01 

6.488E+01 

4, 820E+01 

4.040E+01 

.4500 

2 320 5. 6 i 

• 76 56 

34.9329 

64.9054 

48.2755 

40.6261 

4.862E-05 

4.774E+C3 

3.498E+01 

6.491E+01 

4.828E+01 

4.063E+01 

.460? 

23246.60 

. 76d5 

34.9893 

64.9195 

48.3018 

40.9436 

4.888F-05 

4. 782E+C3 

3.499E+C1 

6.492E+01 

4. 830E+01 

4.094E+01 

.4700 

20249. J 5 

• 766 b 

34.9394 

64. 9496 

48.3299 

40.9891 

4. 892E-05 

4. 784E+C3 

3.499E + 01 

6.495E+01 

4.833E+01 

4.099E+01 

.4800 

23226. 74 

• 7663 

34.9904 

64.9539 

48.4048 

41.0153 

4.904F-C5 

4. 788E+C3 

3.499E+01 

6.495E+01 

4. 840E+0 1 

4.102E+01 

.4900 

23261.73 

.7676 

34.9920 

64.9805 

48. 5162 

41.2520 

4. 916E-05 

4, 792E+03 

3. 4 99 E+Ol 

6.498E+01 

4.852E+01 

4.125E+01 

.5000 

23262. 16 

. 9699 

35.CC51 

65.0048 

48.7771 

41.3566 

4. 920E-C5 

4.794E+03 

3. 501E+C1 

6.500E+01 

4. 878E+01 

4.136E+01 



TABLE III.- SAMPLE OUTPUT - Concluded 


.5100 

20325. 4? 

.7682 

3 5.0111 

65.0207 

48. 9386 

41.3709 

4. 930E-05 

4.797E+C3 

3.501E+01 

6.502E+01 

4. 894E+01 

4.137E+01 

. 52C? 

23342. 60 

.7632 

35. 023° 

65.0360 

49.4466 

41.4382 

4.931E-05 

4 , 8C2E+03 

3.502F+01 

6. 504E+01 

4.945E+01 

4.144E+01 

.5300 

23 3°2« 36 

• 76 37 

35.0285 

65. 0888 

49.4754 

41.6255 

4. 937E-05 

4.8C3E+C3 

3 . 503 E + C 1 

6.509E+01 

4.948E+C 1 

4.163E+01 

. 5 ACC 

2 34+3. 42 

. 7693 

° 5. 0309 

65.C887 

49.4781 

4!. 7*98 

4. 938E-05 

4.8C6E+C3 

3. 503E+01 

6.509E+01 

4.948E+01 

4.175E+01 

.5500 

20453. 1= 

. 7693 

35.0397 

65. 09C1 

49.65C2 

41.7630 

4. 939E-C5 

4.817E+03 

?. 504E+01 

6.509E+01 

4.965E+01 

4.176E+01 

. 56CC 

20567.06 

.7694 

35.0663 

65.0940 

49.9076 

41.9083 

4. 949E-05 

4. 822E+03 

3. 507E+01 

6. 509E+01 

4.991E+01 

4.191E+01 

. 570C 

23505.51 

. 7700 

35.0767 

6 5. ice4 

49.9366 

42.2939 

A.965E-05 

*+. 827 E + C3 

3. 5C8E+01 

6.511E+01 

4, 994E+01 

4.229E+01 

. 5 8CC 

206:4. a: 

.33J-J 

35.1066 

65.1953 

49.9733 

42.4465 

4.967E-05 

4. 832E+C3 

3. 511E+01 

6. 520E+01 

4.997E+01 

4.245E+01 

.59CC 

20627. 54 

.3707 

35.1094 

65. 2085 

50.1443 

42.4562 

4.971E-05 

4.8A6E+C3 

3.51’E + Ol 

6.521E+01 

5.014E+C1 

4.246E+01 

. fcOGC 

206-4.70 

.7708 

35. 1222 

65.2173 

50.2390 

42.5089 

4.992E-05 

4.858E + C3 

3. 512E+01 

6.522E+01 

5.024E+01 

4.251E+01 

.6100 

207 28. 32 

.7721 

35. 1353 

65. 2189 

50.3370 

42.5362 

5.006E-05 

4.859E+C3 

3. 514E+01 

6.522E+01 

5.034E+01 

4.254E+01 

. 6200 

20771.64 

.7726 

35.1549 

65.2355 

50.4117 

42.5571 

5. 009E-05 

4. 865E+03 

3. 515E+01 

6.524E+01 

5.041E+01 

4.256E+01 

.6300 

2386 7. C 7 

.7728 

35.1762 

65. 3074 

50.4586 

42.6173 

5.009E-05 

4. 865 E + C 3 

3.518E+01 

6. 5 31E + 01 

5.046E+C 1 

4.262E+01 

• 6A00 

2390 2.99 

.7730 

35. 1767 

65. 32e5 

50. 4975 

42.7145 

5. 017E-C5 

4. 873E+C3 

3.55 8E+01 

6. 533E+01 

5.050E+0 1 

4.271E+01 

. 6500 

20 9 2 3. 34 

. 7731 

35. 1778 

65.3739 

50.8142 

42.8295 

5.021E-05 

4.875E+C3 

3. 5 18E+01 

6. 537E+01 

5.081E+01 

4.283E+01 

.6600 

20946.96 

.7731 

35.1353 

65. 3987 

50. 8b79 

42. 8 5 1C 

5.035E-05 

4. 977E+C3 

3.519E+01 

6.540E+01 

5.087E+01 

4.285E+01 

. 67C0 

20559.46 

.3734 

35. 207° 

65. 41 C 3 

50.9720 

43.0323 

5.049E-05 

4.9C7E+03 

3.521E+01 

6.541E+01 

5.097E+01 

4.303E+01 

.6800 

23973.91 

.773b 

35. 2281 

65.4253 

51.C506 

43.0889 

5.064E-C5 

4.9C9E+03 

3.523 E + 01 

6.543E+01 

5.105E+01 

4. 309E+01 

.6900 

2 10 7 6. 4 9 

.7739 

35.2440 

65.4308 

51.14C5 

43.2585 

5.070E-05 

4.917E + 03 

3. 524E+01 

6.543E+01 

5.114E+01 

4.330E+01 

.7000 

212:9.02 

.7738 

35. 2489 

65.466b 

51.1463 

43.2990 

5.091E-05 

4.919E+03 

3, 525E+01 

6.547E+01 

5.115E+01 

4. 330E+01 

.7100 

214! 0.C4 

.3740 

35. 2590 

65.5145 

51.3315 

43.5703 

5.093E-05 

4.936E+03 

3, 526E+01 

6. 551E+01 

5.133E+C1 

4.357E+01 

.7200 

21441.65 

.7747 

35. 2694 

65.5305 

51.3512 

43.6385 

5. 102E-05 

4.941E+C3 

3.527E+01 

6. 553E+01 

5.135E+01 

4.361E+01 

.7300 

21557.33 

.7747 

35. 2887 

65. 5587 

51.4224 

43.6685 

5. 104E-05 

4.952E+C3 

3. 529E + 01 

6.556E+01 

5. 142E+01 

4.367E+01 

. 7400 

21538.27 

.7757 

35.3132 

65. 5637 

51.9994 

43.7537 

5.133E-05 

4.954E+03 

3.531E+C1 

6.556E+01 

5. 200E+01 

4.375E+01 

. 7500 

216C3.57 

.7761 

35.3242 

65. 5813 

52. 1941 

43. 826? 

5. 134E-05 

4.957E + C3 

3. 532 E + 01 

6.5 58E+01 

5.219E+01 

4.383E+01 

.7600 

21668.47 

.7771 

35. 3276 

65.5941 

52.2220 

43.9649 

5.137E-05 

4.960E+P3 

3. 533E+C1 

fc. 5 59E + 01 

5.222E+C1 

4.356E+01 

.7700 

21642.47 

.7774 

35. 3331 

65.6293 

52. 3561 

44,0916 

5. 143E-05 

4.965E+0 3 

3. 533E+01 

6.563E+01 

5.236E+01 

4.409E+01 

.7800 

21647.39 

. 7792 

35. 3463 

65.6567 

52.4686 

44.3143 

5. 157E-05 

4.982E+C3 

3.535E+01 

6. 566E+01 

5. 247E + 0 1 

4.431E+01 

.7900 

2 1 7 J C. 64 

.7796 

35.3472 

65.6690 

52.6924 

44. 8445 

5.187E-05 

4.983E+C3 

3. 535E+01 

6, 567E+01 

5.269E+01 

4.4B4E+01 

.8000 

21753.42 

.7799 

35. 3935 

65.8725 

52.9177 

44.8562 

5.213E-05 

4. 9865+03 

3, 539E+01 

6.587E+01 

5.292E+01 

4.486E+01 

.8100 

21841.73 

.7815 

35.4078 

65. 5752 

53.6217 

44.95 62 

5. 227E-05 

4.987E+03 

3, 541 E+01 

6.59SE+01 

5.362E+01 

4.496E+01 

. 8200 

2 1 379. 81 

.7823 

35.4145 

65.9876 

53.9690 

45.0106 

5.229E-05 

4.990E+C3 

3. 541 E+01 

6. 599E+01 

5. 397E+01 

4.501E+01 

.8300 

21920.6 7 

. 7e29 

35.4223 

65. 5889 

54.3090 

45.0521 

5.231E-05 

4.990E+03 

3 . 542 E + 01 

6. 5 99 E+01 

5.431E+01 

4.505E+01 

.8400 

21980.29 

.7829 

35.4377 

66.0218 

54.7639 

45.2284 

5.241E-05 

5.014E+03 

3. 544E+01 

6. 602E+01 

5.476E+01 

4.523E+01 

. 8500 

22044. 85 

.7834 

35. 4553 

66. 0397 

54. 8574 

45.2520 

5. 249E-05 

5. 0 1 4E+03 

3.546E+01 

6.604E+01 

5.4B6E+01 

4.525E+01 

. 8600 

2 204 9. 5 2 

.7840 

35.4576 

66.1011 

55.0630 

45.5357 

5. 3t5E-05 

5. 015E+03 

3. 546E+01 

6. 610E+01 

5.506E+01 

4.554E+01 

.8730 

22105.82 

.7840 

35.4584 

66. 1259 

55.0777 

46.0531 

5.315E-05 

5.052E+03 

3. 546E+01 

6. 613E+01 

5.508E+01 

4. 605E+01 

. 890C 

22149.74 

.7842 

35.4596 

66.1636 

55.1636 

46.3 540 

5.315E-05 

5.C61E+C3 

3. 546E+C1 

6.616E+01 

5.516E+01 

4.635E+01 

. 8900 

22361. 64 

.7866 

35.4680 

66.2794 

55.2793 

46.4542 

5.322E-05 

5.091E+C3 

3. 547E+01 

6.628E+01 

5.538E+01 

4. 645E+01 

.9000 

22399.09 

.7879 

35. 5066 

66.3287 

55,9927 

47.3672 

5. 346E-05 

5. 105E+C3 

3. 551E+01 

6. 633E+01 

5.599E+01 

4.737E+01 

.9100 

22501.05 

.7879 

35. 5382 

66.3354 

56.0377 

47.3681 

5.360E-05 

5. 113E+C3 

3. 554E+01 

6.634E+01 

5.604E+01 

4.737E+01 

.9200 

226?C.O° 

.7900 

35.5440 

66.3553 

56.4434 

47. 6594 

5. 361E-05 

5. 1 1 6E + 03 

3. 554E+01 

6.636E+01 

5.644E+01 

4.766E+01 

.9300 

22732.44 

.7836 

35.5610 

66.3906 

56.5616 

49.2 510 

5. 374E-05 

5. 168E+C3 

3. 556E+01 

6.639E+01 

5.656E+01 

4.925E+01 

.9400 

22767. 4C 

. 7568 

35.5700 

66.4176 

56.6055 

49, 6 6 66 

5.379E-05 

5.227E+03 

3.557E+01 

6.642E+C1 

5.661E+D1 

4.967E+01 

.9500 

23463. 86 

. 79 84 

35.6170 

66.4886 

57,1089 

49. 9788 

5.384E-05 

5.244E+C? 

3.562E + 01 

6.649E+01 

5.711E+01 

4.998E+01 

.9600 

24572.03 

. 8C45 

35.6237 

66. 5745 

57.6279 

50.1261 

5.417E-05 

5.252E+C3 

3. 562 E + 01 

6. 657E+01 

5.763E+01 

5.013E+01 

.9700 

24956.00 

.8094 

35.6383 

66.9452 

57.9681 

50.4681 

5.499E-05 

5.255E + C3 

3. 564E+01 

6.695E+01 

5.797E+01 

5.047E+01 

.9800 

24947.24 

.8111 

35.6588 

66. 9858 

58.9622 

50.7445 

5.512E-05 

5.289E+03 

3.566E+01 

6.699E+01 

5.896E+01 

5.074E+01 

. 9900 

25177.63 

. 83 70 

35.7272 

67. 5607 

59.1320 

51.6823 

5. 549E-05 

5. 379E + C3 

3. 573E+01 

6.756E+01 

5.913E+01 

5.168E+01 

1.0000 

26466. 55 

. 8201 

35. 7657 

67.6549 

60.4932 

51.9226 

5.568E-05 

5. 399E+03 

3. 577E+01 

6. 765E+01 

6.049E+01 

5.192E+01 


TABLE IV.- HISTOGRAM SYMBOLS 


Symbol 

Interval value 

2 

4 

8 

+ - - + 

Add 0 

Add 0 

Add 0 

+ • • + 


Add 1 

Add 1, 2, or 3 

+ + + + 

Add 1 

Add 2 

Add 4 

+ 11 + 


Add 3 

Add 5, 6, or 7 
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